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THE RELATIVISTIC EQUATIONS OF STELLAR STRUCTURE AND EVOLUTION 

KIP S. THORNE 
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California Institute of Technology, Pasadena, California 91125 

ABSTRACT 

The general relativistic equations of stellar structure and 
evolution are reformulated in a notation which makes easy contact 
with Newtonian theory. Also, a general relativistic version of 
the mixing- length formalism for convection is presented. Finally, 
it is argued that in previous work on spherical systems general 
relativity theorists have identified the wrong quantity as "total 
mass-energy inside radius r." 
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I. INTRODUCTION 


The general relativistic equations of stellar structure for zero- 
temperature stars (neutron stars) were first presented in their modern form 
by Oppenheimer and Volkoff (1939). Two decades later the discovery of 
discrete, galactic X-ray sources (Giacconi et al . 1962, Gursky et al. 1963) 
motivated theoretical studies of hot, relativistic neutron stars by Chiu 
and Salpeter (1964), Morton (1964), and Tsuruta (1964); and the huge energy 
requirements of strong radio sources motivated Hoyle and Fowler ( 1963a, b) 
to develop the theory of hot, efficiently convective supermassive stars in 
which, it was soon realized, relativistic effects can be important (Feynman 
1964, Chandrasekhar 1964, Fowler 1964). In response to these developments, 
and others, Bardeen (1965), Misner and Sharp (1965), and Lindquist (1966) 
developed the theory of diffusive heat transfer in relativistic stars, and 
Bondi (1964), Chandrasekhar (1965) and Thorne (1966a) elucidated the rela- 
tivistic version of the Schwarzschild criterion for convection. All of these 
pieces of relativistic stellar theory were put together and combined with, 
relativistic equations for nuclear energy generation by Hameen-Anttila and 
Anttila (1966) and by Thorne (1966b, 1967) to give the currently standard 
version of the relativistic equations of stellar structure and evolution. 

Recently Anna 2ytkow and I began analyzing the structure of red super- 
giant stars with degenerate neutron cores (see the following paper and 
references cited therein). For this purpose the standard relativistic 
stellar equations are unsatisfactory in two ways: (i) they do not make easy 

contact with the standard Newtonian equations; and (ii) they do not include 
a mixing- length formalism for convective energy transport. The purpose of 
this paper is to remedy these defects by (i) translating the relativistic 



equations into a new notation, and (ii) presenting a straightforward rela- 
tivistic generalization of the standard Newtonian mixing- length equations. 

No detailed derivations will be given because the translation from the old 
notation (Thorne 1966b, 1967) to the new is straightforward; and the deriva- 
tion of the relativistic mixing-length theory is identical to the Newtonian 
derivation, if one works in the proper reference frame of a static relativistic 
observer. 

Throughout the paper c.g.s. units will be used; the speed of light c 
and Newton's gravitation constant G will not be set equal to unity. 

IT. FUNDAMENTAL VARIABLES 


As our independent thermodynamic variables we choose the following — 
all of which are determined by measurements using standard, physical rods 
and clocks, in the mean local rest frame of the baryons. After the symbol 
'for each quantity, we indicate its units in brackets. 


3 , 

■ cm -) 


/, . . ( mass of one hydrogen 

(density of "rest mass") = . . ■ / ® ^ 

l atom in its ground state 


X 


/number density) 
of baryons J ’ 


T[°K] 


X. 


= (temperature); 


fractional abundance of nuclear species 1, \ 
by rest mass or equivalently by baryon number/ * 


(la) 

(lb) 

(lc) 


As our independent radial and time variables we choose 


M r [g] 


"rest mass" inside 
"radius r" 


X 


/ mass of one hydrogen 
^atom in its ground state 


/total number of baryons) 

( inside radius r / , 


(2a) 


2 



t[sec] = 


"Schwarzschild Lime) 
coordinate" 1 


\ I time coordinate such that [c>/dt] r is \ 
j = the time- translation Killing vector and . (2b) 
' t is proper time at radial infinity , 


The gravitational field is characterized by three fundamental variables 

which are functions of M and t; 

r 

r[cm] = ("radius") = (l/2or) X (circumference around center of star); 


"total mass inside radius r" — including contributions from rest) \ 
mass, nuclear binding energy, internal energy, and gravity /*' J 


$ r- ff-1 = ("gravitational potential") = 
Lsec J 


* c2 * n [it 


Energy transport through the star is characterized by three quantities, 

/ 

each of which is determined by measurements using standard, physical rods and 
clocks in the mean local rest frame of the baryons at radius r: 


'ergl = ("i ocal luminosity") = /non-neutrino energy being transported across 
sec ~ \ the sphere at radius r, per unit time 


~erg~ l " 
sec “If 


"neutrino luminosity ' 
from nuclear burning" 


same as above, but for neutrino energy) 
produced in thermonuclear reaction ; 
cycles which change the abundances X. i 


erg _ 
sec = ^ 


- \ /same as above, but for neutrino\ 

neutrino luminosity not\ / , ’ , , \ n N 

from nuclear burning" M energy produced by processes . ‘(4c) 

’ 1 which, m time- average, do not 

\ change the abundances X. / 


The complete stellar structure and evolution are characterized by the 
functions p(M r ,t), T(M r ,t), X^M^t), r(M r ,t), ^(M^t), ®(M r ,t), L r (M r ,t), 
b" V (M r ,t), and L° V (M r ,t). 
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III. AUXILIARY VARIABLES 


a) Thermodynamic, Nuclear Burning, and Opacity Variables 


The following auxiliary variables are algebraic functions of the funda- 
mental variables; and like the fundamental variables they are determined by 
measurements using standard, physical rods and clocks in the mean local rest 
frame of the baryons 

P(p,T,X^) te-gil == (total pressure); (5a) 

L cm J 



binding energy of nuclei, per unit\ 
rest mass, relative to hydrogen J 


= 1-2 


m. X. 1 

l l 


i"H A i 


(5b) 


where nu is the mass of atomic species i in its ground state, m^ is the mass 
of atomic hydrogen, A^ is the number of baryons in atomic species i, and c is 
the speed of light; 


n( P ,T,x.) 



= ("specific internal energy") 


| total mass-energy of a sample of ] 
ystellar material, in energy units/ + g(x ) - * 

(total rest mass of the sample) ' i ’ 


(5c) 


p„(p,T,X.) PU = f denSity ° £ total non-gravitati 0 nal\ p(1 _ B / C 2 + It/ C 2 ); (Sd) 
t ’ i' [_ cm 3J y mass-energy, in mass units J ’ 


(p,T,X^) = (opacity) = (Rosseland mean opacity); 


(5e) 



e (p,T,X.) 

nuc' K ' 7 i 7 

e (p,T,X.) 
nv 7 7 i 7 

e (p,T,X.) 
ov 7 7 1 / 

^(P,T,X.) 



/rate, per unit rest mass, at which nuclear \ > 
l burning creates non-neutrino energy ) 7 


/rate, per unit rest mass, at which nuclear^ c 
y burning creates neutrino energy J 7 


I rate , per unit rest mass, at which non-nuclear-\ 
burning processes [processes with no change 
^ in X^] create neutrino energy J 


rate at which the abundance X^ of species i 
changes due to nuclear burning 


(5f) 


(5g) 


; (5h) 


(5i) 


b) Relativistic Correction Functions 

The above auxiliary variables (except B and p ) are all familiar from 
the Newtonian theory of stellar interiors. In the relativistic theory it is 
useful to introduce the following additional auxiliary variables, each of 
which is dimensionless and is unity in the Newtonian limit 


ft = ("redshift correction factor") 

= exp($/ c 2 ); 

(6a) 

If = ("volume correction factor") 

= (l-2GM tr /c 2 r)" 1 / 2 } 

(6b) 

^ /"gravitational-acceleration\ 
= ( correction factor" J 

M + !+3tr 3 p/ c 2 
tr ' 

r= M * 

r 

(6c) 

g = ("energy correction factor") 

s i + (n-B)/c 2 = p t /p; 

(6d) 

H = ("enthalpy correction factor") 

p 

= i + (n - b + p/p)/c . 

(6e) 


In terms of these variables, the general relativistic metric for spacetime 
inside and around the star is 

j 2 2 2 , 2 2 , 2 2/ ,„2 , 2 _, 2 \ r \ 

ds = 4 c dt + if dr + r (dO + sm 9dcp ). (7) 


o 



c) Mixing-Length Variables 


The Newtonian mixing- length theory of convective energy transport is 
readily generalized to general relativity. One need only introduce the 
local proper reference frame of an observer at rest at radius r , and in that 
reference frame analyze, in a manner identical to Newtonian theory, the 
buoyant forces on convective cells and the heat exchange between convective 
cells and their surroundings. The auxiliary variables that enter into such 
an analysis, patterned after Paczynski's (1969) Newtonian variant, are 


cm 1 /local acceleration\ GM r 


cm 

2 = 
t-sec J 


\ of gravity 


■ Sir j 


(8a) 


H_ [cm] == (pressure scale height) = (p/pg)?/ ^ ; 


(8b) 


i [cm] s (mixing length [normally chosen equal to Hp]); 


(8c) 


w s (optical thickness of one scale height) = k p 


(8d) 


c = (specific heat at constant pressure) = * (8e) 


P,X. p^-'p, 


.X. 


'.[■ 


sec 

cm 


G p ^ 1 + co^/3 

s (coefficient of heat exchange) = - - - „ ■ L — 


(8f) 


8crT u 


where a = ac/4 is the Stephan-Boltzmann constant} and 


-I 


( S £n p 

Tin T 


p ’ x i 


(8g) 


In terms of these auxiliary variables, the basic algebraic equations of the 
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mixing- length theory are these: (i) An equation which defines the "radiative 

gradient" 

^ _ /value that (d £n T/ciMj.) t ( ci in p/cSMj. ) t * = dinT/d fnP\ 

rac * \ would have if the material were non- convective ) ’ 

the equation for V ^ follows from equations (3.11-3) and (3. 11- 7a) of Thorne 
(1966, 1967) by straightforward change of notation: 


V rad eb'n 


k L r P 


GM ctT 
r 


?f K& 


L_ + fi . 

tv u j 


Oa) 


(ii) The usual equation for the "adiabatic gradient" 

_ / d lnT \ 

- \btnv) t 


r 2 " 1 


rad 


(9b) 


entropy,X^ 


where is the adiabatic index of the second kind, (iii) A set of four 
coupled algebraic equations which determine the energy flux carried by con- 
vection F conv > the mean velocity of a convective cell ("turbulent velocity") 
v , the gradient associated with a convective cell V 1 , and the actual gradient 
averaged over all convective cells and over the medium through which they 
move, v: 


F 


l6crT 


conv 3npHp ' rad 




F = i C_ pTv fi /hJ (v - V'), 

conv a P K t r P' v ' ’ 

v t 2 = |-s WV Q(v - v,) ; 


(10a) 


(10b) 


(10c) 


(V- V')/(V' - V ad ) = 7 0 v fc . 


(lOd) 


Equations (lOb)-(lOd) have identically the same form as in Newtonian theory 
because their derivation in the proper reference frame of a static observer 
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is identical to that of Newtonian theory. Equation (10a) also has standard 

Newtonian form. It follows from equation (9a) with rewritten as 
2 

(^ conv + F ra d) > an ^ ^ rom analogous equation for the actual gradient 
V in terms of the radiative flux F , 


5 _ *+* A r ad> F t 

GM<rT 4 *** 



Because equations (10) all have the same form as in Newtonian theory, one 

can use the standard technique [eqs. (22)-(27) of Paczynski (1969)] to solve 

them for the four unknowns F , v^, V' and V. 

comr t* ’ 


IV. DIFFERENTIAL EQUATIONS OF STELLAR STRUCTURE 

There are 8+N (where N is the number of nuclear species) differential 

equations of stellar structure for the 8+N fundamental variables p, T, 

r, M. , L , L nV , and L° V as functions of M and t. In these differential 
y tr> ’ r’ r 9 r r 

equations ci/ciM^ acts at fixed t, and d/dt acts at fixed M^. Each equation is 
a translation of the indicated combination of equations from Thorne (1966b, 

1967). 

The equation for as a proper volume integral of p$ translation of 
equation (3.11-1): 

dr/dM r = (lla) 

The equation for total mass-energy inside radius rj translation of 
equations (3.11-2) and (3.11-1): 

c>M tr /dM r = e/r. (lib) 


r- , ' 
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The source equation for the gravitational potential $; translation of 
equations (3.11-4) and (3.11-1): 

r 4rtr p 

The equation of energy generation; translation of equations (3.11-5), 
(3.11-6), and (3.11-1): 


x S(L r R 2 ) 


E nuc ” e ov 


i bn jp_ i Sp 

ft + 2 ft bt * 


The equation for neutrino losses due to nuclear burning; translation 
of equations (3.11-6) and (3.11-1), specialized to nuclear-burning neutrinos 

x 5 (l"V) 


The equation for non- nuclear-burning neutrino losses; translation of 
equations (3.11-6) and 3.11-1), specialized to non-nuclear-burning neutrinos 

i S(L ov a 2 ) 

1 v r ' (. 


e ov * 


The equation for changes of nuclear abundances due to nuclear burning; 
translation of equation (3.8) 

r -1 bx./^t - a. , ( 

i i 

The equation of energy transport; follows directly from the definition 
of V; translation of the mixing- length-generalized version of equations 
(3.11-7) 

S In T/SM = V . 3 InP/dll if £ V ,, 


d In T/SM r = V d In P/dM r 


if 

^rad ~ ^ad> 

if 

^rad > V ad‘ 
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The Oppenheimer-Volkoff equation of hydrostatic equilibrium; translation 
of equations (3.11-3) and (3.11-1) 

lr = -5 w - (1 

This equation must be combined with the equation of state P(p,T,X^) and with 
equations (llg,h) for St/^M^ and c&^/ciM^ to yield cip/dM^. 


V. BOUNDARY CONDITIONS 

Corresponding to each different derivative with respect to in the 
equations of stellar structure there is a radial boundary condition. The 
obvious boundary conditions at the star's center are 

r = M. = L = L nV = L OV = 0 at M = 0 (12a) 

tr r r r r s ' 

(translation of [3.38a]). 

We shall denote the surface values of rest mass, total mass, radius, and 
the total luminosities by 

M = M , M. = M , R = r, L = L , L nV = L nV , L° V = L° V at surface. (13) 

r* t tr* ’ r ’ r ’ r s ' 

At the surface the star's spacetime geometry (7) must match onto the external 
Schwarzs child geometry 

ds 2 = -(1 - 2GM t /c 2 r)c 2 dt 2 + (l - 2GM t /e 2 r)“ 1 dr 2 + r 2 (d© 2 + sin 2 9 dep 2 ) . (l4) 

Smoothness of the match ("continuity of intrinsic geometry of surface") 
requires that $ satisfy the surface boundary condition 

$ = \ C 2 (l - 2GM t /c 2 R) at M = M (surface of star). (12b) 
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Note that the luminosities as measured far from the star — which we denote 

nv ov nv 

=c > £ ) and — are not the same as the surface luminosities L, L , and 

L° V . Rather, they are the surface luminosities corrected for gravitational 

redshif t 

s£/L = £ V /L nV = ^ OV /L° V = (1 - 2GM t /c 2 R) . (15) 

In addition to the boundary conditions (l2a,b) one must also impose 
surface boundary conditions on pressure P and temperature T. If moderate 
errors near the surface are allowable, one can impose the "zero boundary 
conditions" 

p = T = 0 at M r =M (12c) 

(translation of eqs. [3.38c,d]). If higher accuracy is desired one can 
impose the boundary conditions of the relativistic version of the Eddington 
approximation 

L = 4jtR 2 crT\ kP = f- (GM t /R 2 )r at = M (12c') 

(translation of eqs. [3.38c r ,d ']) . For still higher accuracy one can join 
onto a model stellar atmosphere. If the atmosphere is thin compared to the 
stellar radius R, then it can be constructed in the standard Newtonian manner 
using a surface gravity of 

g g = (GM fc /r 2 ) If at M r = M, (16) 

a surface luminosity equal to L, and radial and time coordinates t and t 
related to r and t by 

r =(r-R )r, t"=tft; A = t' 1 = ( 1 - 2GM t /c 2 R ) l ^ 2 . (17) 
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If the atmosphere is not thin compared to R, one can construct it using the 
formalism of general relativistic radiative transfer theory, which is reviewed 
in §2.6 of Novikov and Thorne (1973). In the case of a thin atmosphere all 
spectral features as observed by a distant observer are redshifted relative 
to their rest wave lengths by 

AA/A = (l-2GM t /c 2 R)“2-l. (18) 

VI. SOME USEFUL RELATIONS 


In this section we list several useful relations among the stellar- 
interior variables. 

The sum of the fractional abundances must be unity at all times; 
and consequently, the sum of their rates of change must vanish: 


E X. = 1, 
i 


E a. = o. 

i 1 


(19) 


The total rate of energy release by nuclear burning must equal the rate 
of change of nuclear binding energy 

2 


nuc + £ nv R cit 


- s. _ m. c 

1 dB _ i 

= - E : — a. 


i Vi 


( 20 ) 


see equations (5b) and (llg). 

Th,e rate of change of the total mass-energy inside radius r, as measured 
by an observer there, must be equal to the rate at which matter carries mass- 
energy inward minus the rate at which luminosity carries it outward; 



tr 






( 21 ) 
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This mass-energy conservation law requires some discussion: (i) The time 

derivatives here are taken at fixed radius r, whereas all previous time 
derivatives were taken at fixed rest mass M j the two types of time deriva- 
tives are related by 

(d/c3t) r = + (dM r /dt) r (V^M r ) t • ( 22 ) 

r 


(ii) The operator ft - ^(d/dt) is derivative with respect to the proper time 
of an observer who sits at rest at radius r; see equation (7). (iii) 
ft" 1 (£M r /dt) r is the locally measured rate at which rest mass flows inward 
across radius rj and (dM, i^/Bt) ^ is the rate of inflow of rest mass 

plus enthalpy in mass units. Enthalpy appears in the' conservation law 
rather than energy rather than e) for the same reason as it appears in 

the Bernoulli equation: in moving matter, pressure (the difference between 

pW and pg) transports energy 

(energy flux) = (pressure) X (velocity). 

(iv) <l/c 2 )(L + l“ v + L° V ) is the locally measured rate at which mass-energy is 
transported outward by neutrinos, photons, and diffusive heat flow. (v) Since 
M tr is the total mass-energy inside radius r, one would have expected the 
left side of equation (21) to read (dM /dt') r — i.e., one would have 
expected the y to be absent. The presence of y suggests to me that 
relativity theorists such as Misner, Thorne, and Wheeler (1975) should not 
have given the name "total mass-energy inside radius r" to M tr . Rather, 
the quantity 


tr = ( c2 / G ) r = (c 2 /G)r 


1 - (1 - 2GM tr /c 2 r)2 ^ 


(23) 


2 2 

— + b GItf /c r in Newtonian limit 

tr ^ tr 
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should have been identified as total mass-energy inside radius r because 


it satisfies 




-1 


-«- 1 fr)-- 1 © 


V - \ (L r + L^ V + L° V ). ( 21 ') 


Out of deference to established convention I suggest that people retain the 
name "total mass-energy inside radius r" for M , but keep in mind that it 
is a misnomer. 

The equation of mass-energy conservation (21) can be derived from the 
equations of stellar structure by first deriving the relation 


a 

dF 


r i _ /oM \ /oM. \ 1 

r (L. + C V +L° V )|* - (-gf-J va + (s^j Vft = 0, 

' r ' / r 


i-c 


r\ V I - 1 - 1 

2 r r 


(24) 


where d/dr acts at fixed time t, and by then invoking the boundary conditions 


T T nv T ov M .. _ 

L -L = L = M = M ■ = 0 
r r r r tr 


at r = 0, 


A derivation of equation (24) proceeds as follows: (i) By combining equations 

(lla,d,e,f), (20), and (6d) derive the relation 


d 

dF 


4 d. + L« + i»V' 


*- c 


r r 


= -4jtr plrfo 


de\ _ p/p 2 / dp\ 
dt I ' 2 l dt L. 

/M c \ JH 


(ii) Use equations (22) to convert from time derivatives at fixed M r to time 
derivatives at fixed rj and then use equations (lla, c, i) and (6a, d, e) to obtain 


\ (L +L nV + L° V )g 2 
2 r r r 1 


-c 


-4itr pTrft 


ds\ P/p 2 / dp\ 
dE " 2 \dtj 

/ r c ' * r J 


/dM \ 

+ \~df) dF * 


(25a) 
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(iii) Use equation (lla) to derive the relation 


s? (sf) = (”^r) ^ + wVft [ a ^f- ] • ( 25b ) 

* r r -* r 

(iv) Use equations (6a,b, c,d,e) and (lla, b, c) to derive the relation 

(d/dr) (r ft) = 4jt(G/c 2 )rp^r 5 ^ 

and then use equations (lla,b) and (6b) to obtain 

& [(?r ) r «] - ^ H [^] r + r r | • < 25 ‘> 

(v) Finally, combine equations (25a, b,c) and (6d,e) to obtain equation (2^). 
Note that the equation of mass-energy conservation (21), when evaluated far 

outside the star, just says that if the rest mass of the star is held fixed 
then its total mass-energy decreases at a rate given by the photon and neutrino 
mass-energy losses 

dM t /dt = - (l/c 2 )U + =£ nV + =£° V )* (26) 

VII. SUMMARY 

Coordinates M^, t for the stellar interior are defined in equations (2a,b). 
The star’s structure is described by 8+N (where N is the number of nuclear 
species) fundamental variables p, T, X., r, M fcr , $, L r , l£ V , and L° V , which 
are functions of M r and t, and which are defined in equations (l), (3), and 
(^). These 8 + N variables satisfy the 8 + N differential equations of stellar 
structure ( 11), subject to the radial boundary conditions (12). The differen- 
tial equations (ll) contain a number of auxiliary variables, which are 
algebraic functions of the fundamental variables, and which are defined in 

15 



equations (5)-(l0). Quantities which characterize the surface of the star, 
its external gravitational field, and the radiation which leaves the star 
are described by equations (13)-(18). Several useful relations among the 
stellar variables are given in equations (l9)-(2l) and (26). 

This version of the equations of stellar structure and evolution reduces 
to the standard Newtonian version when one sets the following relativistic 
correction factors to unity: R, If, &, 2, % in the interior; (l~2GM t /c R) 

and ( 1 - 2GM t /c 2 r) at the surface and in the exterior. 
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ABSTRACT 

Stars with massive envelopes (M env ^ 1 %) and degenerate neutron 
cores (M core — 1 %), R core ^ 10 km) are analyzed theoretically: General 
relativistic equations of structure are derived under the assumptions 
of hydrostatic and thermal equilibrium, spherical symmetry, no rotation, 
and no magnetic field. Numerical models are constructed, and analytic 
expressions are derived for the stellar structure in various interior 
regions. It is argued that all nonrotating, equilibrium models prob- 
ably resemble qualitatively those constructed in this paper. Brief 
discussions are given of the stability and evolution of the models, and 
of prospects for identifying such stars observationally. 

Viewed externally, our models are extreme M supergiants (l~j3x 10^ 
to 1.3X105 Io, T photo here ~ 2600 to 3100 K, R photo „ p u ere - 1000 %,) . 
The large, diffuse envelope of each model is separated from its compact 
core by a thin (~ ^0 meter) energy-generation layer called the "halo." 
The envelope convects from the outer edge of the halo all the way out 
to the photosphere. Matter contracts from the envelope through the 
halo and into the core at a rate of ~ 1x10“® P^/yr. The contracting 
matter releases its gravitational energy and burns its hydrogen and 
helium while passing through the halo. When the envelope mass exceeds 
~ 10 1^), the hydrogen-burning shell occurs at the halo-envelope inter- 
face, and the products of hot (T 1 X 10® K) nonequilibrium hydrogen 
burning are convected directly from the burning shell out to the photo- 
sphere, where they should be observable. 
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I. INTRODUCTION AND OVERVIEW 


a) Stars with Neutron Cores Compared with Stars with White-Dwarf Cores 


This is the first of several papers devoted to the question "What are 
the possible equilibrium states for a star consisting of a massive nondegen 
erate envelope surrounding a degenerate neutron core?" 

The analogous question, "What are the equilibrium states for a star 
with a massive, nondegenerate envelope surrounding a degenerate-electron 
(white-dwarf) core?" has a well-known answer: Such stars are red giants 

which reside near the Hayashi track of the H-R diagram. In these stars mat 
ter continually, but slowly, flows from the inner regions of the envelope 
onto the outer regions of the core, passing through one or more nuclear 
burning shells as it flows. The inflow releases nuclear and gravitational 
energy, converting it into stellar luminosity L : 


L = L + L , L 

nuc grav nuc 


' 2 

Me Q, 


L 

grav 




3 


(1.1a) 


Q - 0.007 


GM /R c 2 
c c 



(1.1b) 


(Here M is the rate of mass flow into the core, Q is the efficiency of 

2 

nuclear burning for converting rest mass into thermal energy, and GM^/R^c 
is the analogous efficiency of gravitational contraction with and R^ 

the core mass and radius.) 

For the case of a star with neutron core, one might expect a similar 
answer: Red giant star near the Hayashi track; gradual inflow of matter 

from envelope to core; formula (1.1a) for energy generation again valid, 
but now with 


1 




Q - 0.007 , 


(1.1b’) 


GM /R c 2 - 0.15 . 

c c 

The enormous strengthening of the gravitational potential, GM^/R^ , 
when the white-dwarf core is replaced by a neutron core, has two conse- 
quences: (i) The relative roles of nuclear burning and gravitation as 

sources of luminosity are reversed: 


L /L - 0.99 , L /L - 0.01 for white-dwarf core , (1.2a) 

nuc grav 


L /L - 0.04 , L /L - 0.96 for neutron core . (1.2b) 

nuc grav 


(ii) The timescale for marked evolution of the star is much longer 

in the neutron case than in the white-dwarf case, if one compares stars of 

similar luminosities: 


T . 1M /M L /0.007 c L , 

neut „ o neut K w. d . _ 20 w * » 2 . 


(1.3) 


w.d. 1M /M , L /0. 15 c 
© w.d. neut 


neut 


b) Qualitative Overview of the Internal Structure 


The above discussion is corroborated by the detailed stellar models 

that we shall construct in this paper — so long as the total mass of the 

star is < 10 M (we shall call such, stars "giants"). For M > 10 M 
~ 0 0 

("supergiants") our models have convective envelopes that extend all the 

way into the hydrogen-burning shell. As a consequence, most of the burned 

material is recycled back into the envelope, rather than being passed on 

into the core; the relative importance of nuclear and gravitational energy 

generation is reversed back to the white-dwarf-type situation, » 

L ; the evolution of the star is dominated by chemical changes in the 
grav 


r. 
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envelope rather than by growth of the core; and the evolution timescale 
, is comparable to the white-dwarf-core case. 

Except for location of the hydrogen-burning shell and its resulting 
influence on the star’s evolution (giant versus supergiant), our stellar 
models all have similar structures. Figure 1 depicts their common struc- 
ture, and defines a number of terms ("envelope", "knee", "halo", "core", 
...) which we shall use throughout this paper in discussing our models. 

The stellar structure depicted in Figure 1 is very peculiar; many 
of its features are unique to stars with neutron cores, and violate in- 
tuition based on studies of more normal stellar models. For example: 

(i) In no other type of stellar model yet constructed does a single con- 
vection zone link the photosphere to a nuclear-burning region, (ii) The 
region between the core and the base of the convective envelope is nearly 
isothermal and has a total thickness of only 'v 40 meters; we call this 
region the star’s "halo". (iii) All of the gravitational energy release 
occurs in the upper regions (< 20 meters) of this halo. (iv) In giant 
models this halo contains both the hydrogen- and helium-burning shell 
sources, each with thickness < 5 meters; in supergiants the hydrogen- 
burning shell overlaps the envelope, so the halo contains only the helium- 
burning shell. 

c) Observable Features of the Models 

Unfortunately all of these extreme halo conditions are thoroughly 
hidden from the prying eyes of the astronomer by the huge, tenuous, red- 
giant envelope. The envelope acts as a buffer: Consider two stellar models 
with the same core mass, envelope mass, and total luminosity, but with 
different cores (white-dwarf versus neutron). Imagine comparing these 
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models by swimming outward from the core through the envelope to the 
photosphere. The differences one would see are enormous near the core; 
but they would gradually die away as one moves outward through the envelope. 
At the photosphere only one tiny difference would remain: the star with 
neutron core would be slightly redder, by A log T ^ << 0.1 . 

Put differently: aside from chemical composition (see below), the 
only distinguishing external feature of our models with neutron cores is 
their extreme redness: because they sit precisely on the edge of the 
Hayashi forbidden region, they must be the reddest stars in the universe; 
but they will be redder than stars at the tip of the normal giant branch 
by only a very slight amount, A log T ^ < 0.1 . This difference is so 
slight that it will get lost in other effects (reddening by circumstellar 
material and interstellar material, differences in chemical composition 
causing differences in T h » uncertainties in values of molecular opacities 
and convective mixing lengths, etc.). Hence, this redness difference is not 
a good "handle" to use in observational searches for stars with neutron cores. 

Thus far our model building has yielded only one good observational 
handle — and we are not yet sure of its details: In our supergiant models 
convection should carry the products of hydrogen burning directly from 
the nuclear -burning shell to the photosphere. The hydrogen will be burned 
by a hot (T - 1 x 10 9 K) , nonequilibrium CNO-Ne reaction network, and pre- 
sumably will produce very peculiar relative abundances of various catalyst 


18. 17. 16. 13 


isotopes (0, 0, 


0 , 


12 


C, etc.). It may be possible to measure thes 


abundances in the photosphere by observational studies of molecular band 
spectra — e.g., rotational bands of carbon monoxide, vanadium oxide and 
titanium oxide. In collaboration with Michael Newman we are now calculat- 
ing the details of the nuclear reaction chains and the resulting abundances; 
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we shall publish them in a subsequent paper in this series , 

It is conceivable that our models may experience instabilities that 
do not occur in white-dwarf-cored stars with massive envelopes — and that 
the effects of these instabilities might be discemable observationally . 
However, we have not yet undertaken detailed stability analyses of our 
models . 

Our preliminary, crude studies of stability suggest that the envelopes 

of our models might be unstable against complete disruption for M < 3M 

@ 

when M — 1 M However, it seems quite possible that our stars are 
core © 

stable against disruption if M > 5 M q , and in this case live for 3 O' 7 to 

TO 8 

10 years . 

Although a red giant of given luminosity may live 20 times longer if it has 
a neutron core than if it has a white-dwarf core, giants with neutron cores 
may well be much less abundant in the universe than giants with white-dwarf 
cores: When massive stars form neutron cores by gravitational collapse, 

their loosely bound, tenuous envelopes probably get ejected. If so, then 
the only way the neutron core can become a red giant is by acquiring a new 
envelope — and the only place this is likely to happen is in a very close 
binary system, by supercritical mass transfer from a companion or by a 
cannibalistic sinking into the companion's center and eating of the compan- 
ion's core. Recently Ostriker and Paczynski (1975) have speculated about 
such events. 

d) Previous Work on Stars with Neutron Cores 

We are aware of the following previous work on stars with neutron cores: 
(i) In the 1930 's a number of people speculated about the structures and 
stellar-evolutionary roles of such stars, but no detailed analyses were 
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carried out and no firm conclusions were reached; see, e.g., Gamow (1937), 
Landau (1937), Oppenheimer and Volkoff (1939). For example, Landau (1937) 
noticed the enormous efficiency, GM^/R^c - 0.15, with which contraction onto 
a neutron core can liberate energy; he proposed that this might be the source 
of the luminosity of the sun and other stars; and he suggested that one try 
to build stellar models of this type. (Presumably nobody tried because 
shortly thereafter nuclear burning was recognized as the true energy source.) 
(ii) Murray Gell-Mann tells us that in the early 1950’s Enrico Fermi specu- 
lated that stars with neutron cores would be red supergiants; however, so 
far as we have been able to learn, Fermi never published anything on this 
subject. (iii) Zel’dovich, Ivanova, and Nadyozhin (1972) studied the con- 

_5 

traction of small-mass envelopes (M <10 M ) onto neutron stars. They 

env - o 

found a neutrino luminosity far greater than the photon luminosity; and 
they speculated that, by analogy, stars with neutron cores and massive en- 
velopes might be unstable against collapse of the envelope onto the core, 
with the collapse energy being carried off by neutrinos. We shall argue 
later (§VI below) that our models do not suffer this "neutrino -runaway in- 
stability." (iv) Stothers and Cheng (1974) speculated that the envelope of 
a star with a neutron core would be- rapidly ejected by a secalar instability. 
Our studies (§VI below) suggest that this might be correct fur low-mass en- 
velopes (M^ < 2M q ), but that more massive envelopes might be -stable against 

disruption. ' (iv) Paczydski (private communication, 1973) suggested that one 
of us (ANZ) try to construct stellar models with neutron cores, and we de- 
cided to collaborate on the project. We published a brief account of our 
results as Thome and Zytkow (1975). (v) Ostriker and Paczynski (1975) 

speculated on the role of such stars in the evolution of close binary systems 
(see above) . 
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In the last section of this paper we shall list a number of further 
investigations that are needed. 

e) Notation Used in This Paper 

We here summarize for future reference those mathematical symbols 
which are used in more than one place in this paper, and we give reference 
to equations which contain further details. Equation numbers beginning 
with T are in the accompanying paper (Thorne 1976) ; those beginning with 
A are in the appendix of this paper. We list first the "main symbols" 
and then the "sub and superscripts". 

MAIN SYMBOLS 


conv 


a 
B 
c 
C 

e 

F 

g 

G 

G(T) 

& 

H 

% 

k 


radiation constant appearing in p ra( j = "3 aT 

nuclear binding energy per unit rest mass; eqs . (T,5b) and (2.23). 
speed of light. 

12 

abundance of carbon ( C) by mass, 
relativistic energy correction factor; eq. (T, 6 d). 
energy flux carried convectively ; eq. (T,10). 
local acceleration of gravity; eq. (T, 8 a). 

Newton's gravitation constant. 

Sampson's Klein-Nishina correction factor for electron-scattering 
opacity; eq. (2.32c). 

relativistic gravitational-acceleration correction factor; eq. 

(T, 6 c) . 

scale height [H p eq. (T, 8 b); Hp eq. (i|-.5)$ eq. (4.20b)]. 

relativistic enthalpy correction factor; eq. (T, 6 e) [U eq. (4.2) ] 

• ■ 8 

Boltzmann constant. 


^^t:ox>uc5Ibhaty of y.j.- 
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1. 


mixing length; eq. (T,8c). 


L V 

r 

M 

M 

t 

M 


total non-neutrino luminosity as measured at photosphere and by 
observers far from star; eq. (4.8b) [ L g rav ec *‘ (4»8a)]. 

non-neutrino luminosity as measured at radius r inside star; 
eqs. (T,4a) and (2.43) [L auC eq. (4.16); L ardt eq. (4.4b); 
l/ ad eq. (4.3)]. 

total neutrino luminosity as measured at photosphere [L° V eq. 


.nV 


(T,13); L eq. (1,13)] 


oV 


total neutrino luminosity as measured at radius r [L eq 


(T, 4c) ; L aV eqs. (1,4b) and (2.42a)]. 
mass of hydrogen atom . 

total rest mass; eq. (T,13) [M r eq . (1,2a)]. 

total mass-energy; eq. (l,13)[M tr eq. (T,3b)]. 

rate of inflow of rest mass from envelope to core; eqs. (2.14) 
and (2.44). 

P pressure; eq. (T,5a). 

r and R radius; equal to (1/2tt) x (circumference); eq. (1,3a) and (1,13) 


B 

t 

T 

v 

X 

X. 

x 

Y 

y 

z 


relativistic redshift correction factor; eq. (T,6a) and (2.36) 
[B c eq. (2.42b); B R eq . (2.42b)]. 

Schwarzschild coordinate time; eq. (T,2b). 

Q 

temperature; eq. (T,lb); 1^ E T/10 K; = kT/1 keV. 

locally measured velocity [inflow velocity y eq. (2.19); 
turbulent velocity v eq. (1,10c)]. 

abundance of hydrogen (^H) by mass. 

abundance of nuclear species i by mass; (1,1c). 

abundance of helium (^He) by mass. 

electron-positron pair parameter; eq. (A. 8) . 

1-X-Y ; abundance of "metals" by mass. 
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a 


a. 

i 


Y, 


Yi 

V 

V 


ad 


rad 


luminosity parameter; eq. (2.10); except in Table 3 where a 
is the ratio of mixing length to pressure scale height, 
a = l Jr . 

t p 

nuclear reaction rate for species i ; eq. (T,5i). 


P /P. 
g 


crit 


1 - L /L 
r i 

3 / (1—3 ) = P Vp .. 

g g g rad 


S L /(1-B L ) = L r 


g 

crit 


/L - 1. 
r 


"actual gradient", d in T/d in P; eq. (T,10). 


adiabatic gradient, eq. (T,9b). 
radiative gradient; eq. (T,9a). 

energy generation rate [e nuc eq. (T,5f); eq. (T,5g); 
e oV eq. (T,5h) ] . 

opacity; eq. (T,5e). 


K__ opacity due to scattering of photons by electrons and positrons 

cS 

eq. (2.32). 

K deg e opacity against heat transport by degenerate electrons; eq. 
(2.31) and (2.8). 

y mean molecular weight; eq. (4.9a) £y eq. (2.30); 

y^ on eq. (2.29a)]; except in Appendix where y is chemical 
potential. 

II specific internal energy; eq. (T,5c). 

6 “"3 

p density of rest mass; eq. (T,la); = p/10 g cm ; 

i n ° 

P 1Q = p/10 g cm . 

P Stef an-Boltzmann constant; 0 = ac/4. 

T optical depth measured from the star's surface inward. 

$ gravitational potential; eqs . (T,3c) and (2.34). 

SUBSCRIPTS AND SUPERSCRIPTS 

c outer edge of core (point where electron degeneracy sets in; 

p - 10 6 g/cm 3 ) . 
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c 

cc 

CNO 

crit 

e 

env 

g 

grav 

h 

H 

He 

i 

in 

ion 

K 

m-i 

nuc 

nv 

ov 


o-m 

P 

pair 

ph 

r 


carbon; or at the center of the carbon-burning shell (point 
where C = 0.5). 

12 _ 12 ^ 

C+ C reaction network. 

CNO reaction network for hydrogen burning, 
critical luminosity, 
ionization electrons, 
envelope of star. 

gas (plasma; everything except radiation) . 

gravitational. 

halo of star. 

hydrogen; or at the center of the hydrogen burning shell (point 
outside which half the nuclear energy release has occurred). 

helium; or at the center of the helium burning shell (point 
where Y = 0.5). 

nuclear species i . 

mass inflow from envelope to core. 

ions. 

knee of star. 

at the interface between the middle and inner regions; eq. (2.2). 

non-neutrino energy from nuclear burning. 

neutrino energy from nculear burning; eq. (T,4b). 

neutrino energy from processes other than nuclear burning net- 
works; eq. (T,4c). 

at the interface between the outer and middle regions; eq. (2.1). 
pressure. 

electron-positron pairs, 
photosphere of star, 
measured at radius r . 
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rad 

s 

t 

3a 

+ 


radiation. 

sound. 

turbulence (convective motion) ; except in and M where 
t means "total". 

3a reaction network for helium burning. 

electrons (including ionization electrons and pair electrons), 
positrons. 


f) Outline of Paper 

In §11 we lay down the physical and mathematical foundations for 
the construction of models of stars with neutron cores. Section III is a 
series of graphs and tables displaying the details of our numerical models. 
In §IV we discuss and analyze analytically our "giant models;" and in §V 
we do the same for our "supergiants." In §§VI and VII we discuss briefly 
the stability and evolution of our models. Finally, §VIII is a list of 
topics which need further investigation. 

II. FOUNDATIONS FOR OUR MODEL BUILDING 

In this section we describe the assumptions, equations and numerical 
techniques that underlie our computer-generated models and underlie our 
analytic approximations to them. 

We begin by demanding that our models be spherically symmetric, 
nonrotating, and devoid of magnetic fields, and that they be in slowly 
evolving equilibrium states (evolution timescale long compared to hydro- 
dynamic and thermal timescales). 

Because of the strength of gravity near the neutron core, we ask 
that our models be general relativistic rather than Newtonian — except that 
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Newtonian analyses suffice in the outer region of the star and in order- 
of-magnitude estimates of effects. 


Partition of Model into Three Regions 


In our computer calculations we divide each model into three regions 

(see Fig. 1). The "outer region" contains the atmosphere, the photosphere, 

and the static part of the envelope — i.e. , that portion of the envelope in 

which mass inflow has negligible effects. The " middle region " contains the 

inflowing part of the envelope, the halo, and the outermost layers of the 

core where the carbon-burning shell is located. The "inner region " is the 

entire core, except its outermost layers. 

The boundary between outer and middle regions, r 0 _ m » occurs where 

the inflow first begins to influence the local luminosity ; this happens 

when the enthalpy II + p/p and/or the gravitational potential GM. /r becomes 

2 

larger than ^0. 003 of its maximum value (^0.1 c ). Thus, we arbitrarily set 

/ n p GM 

r = f that radius at which ~tt H r- + — ~ = 3 X 10 

o-m \ ILL 

\ c pc c r 



(All symbols used here are explained in §I.e.) 

The boundary between the middle and inner regions, r^^, occurs where 
nuclear energy generation is no longer significant. In our models more than 
99 percent of all energy generation is by gravity and by thermonuclear 
hydrogen burning, so it is not very necessary to include the effects of 
helium, carbon, or further nuclear burning stages. However, to see what 
their effects may be, we have included helium burning and carbon burning. 

It turns out that the carbon burning is complete by a density of 
p - 1 X 10 g/cm . Therefore, we choose 
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( 2 . 2 ) 


r 

m-x 


(that radius at which 


p = 3 X 1.0^g/cm^) . 


b) Structure of the Inner Region 


In the inner region the high density enforces degeneracy and thereby 
guarantees that the hydrostatic structure (p,P,r as functions of 

M r ) is decoupled from the thermal structure (L r and T as functions of M ) , 
The massive envelope of the star can influence the hydrostatic struc- 
ture of the inner region in only one way: by its weight, which squeezes the 
inner region to a pressure and density, at given M , that are higher than 
for a bare (envelope-free) neutron star. This compressional effect can be 
evaluated by integrating the equation of hydrostatic equilibrium outward 
through the star (throughout this subsection we use Newtonian theory be- 
cause 50 percent accuracy is adequate) : 


M 


M 


P (M ) = 
r 


dP JM 
-j T7- dM 
dM r 


M 


M 


r GM 
= I — ~v dM . 
> 4irr 4 r 


(2.3) 


The fractional contribution of the nondegenerate envelope and halo (region 

g 

with p < 10 ) to the inner-region pressure is 


(AP) 


env 


JL 

P 


GM dM 


10 


6 4TTr 


4 dp 


dp 


(2.4) 


The amount of envelope and halo matter below r = 2 X lQ^km turns out to 

be <10 "^M (see Tables 1 and 2), which is far less than the amount of core 
~ o 

8 8 

matter between p = 1 x 10 and p = 3 x 10 ; hence, in evaluating expres- 
sion (2.4) we can ignore the envelope matter at r < 2 X lO^km — i.e., we 
can regard the envelope as a mass M < 10 M q residing at r > 2 X 10^ km: 
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(AP) 


env 


<1 

~ P 


G(10 M )' 
© 


1 X 10^dynes/cm^ 


4tt( 2 X 10 4 km) 4 


1 X 10 


m-i 


23 


1 X 10 
1 X 10 


23 

26 


« 1 


(2.5) 


This result allows us to conclude that the envelope has no significant in- 
fluence on the hydrostatic structure of the inner region; the inner region 
will have the same hydrostatic structure as a bare (envelope-free) neutron 
star. 

Turn next to the thermal structure of the inner region. At densities 

above p = 3 X 10 11 ("neutron-drip point") the heat conductivity is so high 

that the star is very nearly isothermal [T — const in Newtonian theory; 

T (ft - const in general relativity]. Almost all of the core mass is contained 

in this isothermal region of the core. Between this isothermal core and the 

halo (3 X 10 11 > p > 10 6 ) is a thin "insulating layer" of degenerate-electron 

matter which thermally isolates the core from the rest of the star. We have 

arbitrarily placed our middle-inner region dividing line r m _ i in the center 

8 

of the insulating layer, at p = 3 X 10 . 

Let us estimate the maximum heat flow that the insulating layer can 

support. For ease of computation we shall confine attention to the region 

7 11 

in which the electrons are fully relativistic, 10 < p < 3 X 10 — i.e., 

6 7 

we shall ignore the outermost part of the insulating layer, 10 < p <10 . 

Our estimate relies on the following equation, which is a combination of the 

(relativistic) equation of diffusive heat transfer (eqs. [2.20h] and [T,9a]) 

14 

and the relativistic— degenerate— electron equation of state P = (4.89 X 10 

, / 2 . , , -3 n 4/3 # 

dynes/cm ) (p/g cm ) : 
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d Jin T 
d £n p 


K 


L 

r c 


4/3 


12tt GcM 


= 108 


fC 


10 


M 


tr 


rad 


cm 2 /g 


c M 


tc 


( 2 . 6 ) 


[Here, because the insulating layer is very thin in radius and mass 

(Ar/r << 1 , AM /M << 1), the temperature redshift effect (1 - S/V term 

in V ,) has been ignored, and M has been set equal to the total mass 
rad tr 

of the core M^; also the approximations U - 1 , » anc ^ 

If - lil 1 ~ 61 ^ have been used.] The energy transport is by electron con- 

c 

duction; and the dominant resistance to the conducting electrons in the 
relevant temperature-density regime 

10 7 < p < 3 X 10 11 , 10 8 < T < 10 10 (2.7) 


is electron-electron scattering above the ion-crystal melting temperature 

1/3 

(T_ > T_ n, 1,8 p- _ ' ), and electron-phonon scattering below the melting 
9 9melt 10 

1 / 3 

temperature (Tg < ^ 1 * ^P^q ' )» see flowers and Itoh (1975). In 

these two regimes the computations of Flowers and Itoh give (see their 

Figures 12 and 6 for electron conductivity, which is related to opacity by 

k . K .....= 4ac T 3 /3p) : 
opacity conductivity 

K = (3.9 X 10" 6 cm 2 /g) T 9 4 p 10 ‘ 2 for T g >1-8 P 1Q 1/3 , (2.8a) 

K = (1.11 X 10 _5 cm 2 /g) T 9 3 p 1 “ 5/3 for T g < 1.8 P 1Q 1/3 . (2.8b) 

When inserted into equation (2 . 6) these opacities give 


dT /d Jin p 


10 


-2/3 1/3 

0.35a TnP-irv in molten region, T_>1.8 p, n , 

9 10 9 10 (2 ^ a) 


-1/3 


dTg/d £n p 1Q = a p 1Q in crystaline region, Tg < 1.8 p 


L M _i 
_ , r . , tc. x „ 

a “ 830 L M ' **c * 

e o 


1/3. 


10 


(2.9b) 

( 2 . 10 ) 
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Equations (2.9) and (2.10) set the scale of allowable heat trans- 
fers L through the insulating layer: T Q cannot change by more than 

r y 

a factor 3 , as one traverses the insulating layer, because of the fol- 
lowing: (i) Core neutrino losses keep the isothermal core (and thence its 

11 

outer boundary, p = 3 X 10 ) at a temperature T < 2 . (A neutron star 

y 

cools by neutrino losses to T„ <6 in 12 hours and to T n < 2 

J 9core 9core 

in one year; see Tables 8, 9, and 10 of Tsuruta and Cameron [1966] and 

Fig. 1 of Tsuruta et al. 1972.) (ii) The outer edge of our insulating 

layer, p = 10 , has T Q - 0.5 to 1.0; see Fig. 2, (iii) Neutrino losses, 

y 

which vary as T™ with n > 9/2 in our insulat ing layer (Beaudet, Petrosian, 

and Salpeter 1967), will hold the temperature below T^ - 3 throughout the 

insulating layer. These constraints on T , together with equations (2.9), 

require |a| <1 nearly everywhere in the degenerate electron surface 

layer — and, in fact, |a| < 1/8 in most places including our middle-inner 

8 

region interface, r = r . and p = 3 X 10 . Hence, 

m-i 

|L | < (100 L ) (M /M ) R _1 at p = 3X10 8 . (2.11) 

1 r 1 ~ © tc. © c 


This heat transfer is negligible compared to the star's total luminosity 
L - 10 5 L . 

Nuclear burning of inflowing matter will generate heat in the insulat- 


ing layer of giant models at a rate 


/ dE\ _ /total luminosity \ / GM tc\ ~^- 

\ dt / “ \6£ star, ^5x10 4 L J X \ c 2 R / 


efficiency, 0.0008 of 
mass-to-energy conversion 
when oxygen bums to iron 


) 


'u 300 L (2.12) 

0 

Electron conductivity cannot carry away much more than' 1 -' of this 
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energy; the rest must be carried ’off by neutrinos. 

g 

The above estimates show that the inner region (p > 3 x 10 ) is 
extremely well decoupled from the middle and outer regions, both hydrostat- 
ically and thermally. Its structure and thermal evolution are essentially 
the same as for an isolated (envelope-free) neutron star — and, thus, they 
are not of interest to us here. Henceforth we shall restrict attention to 
the middle and outer regions; and in calculating their structures we shall 
replace the inner region by the "insulation boundary conditions" 

) at r = r . , (2.13a) 

1 / m-x 

at r = r . . (2.13b) 


... .. , . / values for a "bare . . 

(M ,M , and r) = I ' ^ ^ - . .,.o at r = r 

r tr \neutron star at p=3xl0°/ m-x 


m-x 


c) The Outer Region: Physics and Computational Methods 

The outer region includes the atmosphere, the photosphere, and the 
static envelope; see Fig. 1. Our numerical models for this region were 
generated using Paczynski's (1969) computer program "GOB", which calcu- 
lates static stellar envelopes with extended atmospheres using inward 

-12 3 

integrations that begin, in our case, at a density p = 1 x 10 g/cm . 

Each static envelope constructed by GOB can be characterized by the star’s 
total (non-neutrino) luminosity L and mass , the photospheric temper- 

ature T ^ t ^ ie envelope's nuclear abundances — assumed equal to the 
photospheric abundances X Y ph. 5 ^ph * 

The physics and equations which go into the outer region integrations 
are spelled out by Paczynski (1969) . In brief, the physics is this: 

(i) Newtonian equations of stellar structure with luminosity constant 
throughout and with the standard mixing-length formalism for convection; 
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(ii) a simple gray atmosphere model based on the Eddington approximation 
with corrections to account for the "1/r " dilution of the outgoing radia- 
tion, which can be important in extended atmospheres; (iii) an opacity 
table for composition X= 0.7, Z = 0.03 (Paczynski 1970a), 

which is interpolated from the Cox-Stewart (1968) opacities and augmented 
by an approximation to Auman’s (1967) H^O opacity; (iv) an analytic equa- 
t'ion of state including contributions from H^, H, He, H , He , He , free 
electrons, and radiation. 

d) The Middle Region; Physics and Computational Methods 

In the middle region, which we analyze with care, general rela- 
tivistic effects can be important. Therefore, our numerical computations 
utilized the general relativistic equations of stellar structure — which are 
presented in the preceding paper (Thorne 1976; equations in this paper 
are denoted by a T ; e.g., eq.[T,lla]). 

The middle region acts as a conduit through which mass flows from 

the outer region to the inner region. At any given time the total mass 

-8 8 
in this conduit is ^10 M (cf. Tables 1 and 2), which is 'vlO times less 

o 

than the mass in the reservoirs (outer and inner regions) at its two ends. 
Assuming that the star is stable, this huge mass contrast guarantees that 
the rate (per unit Killing-vector-defined coordinate time t) , at which 
rest mass flows inward across a surface of radius r , is independent of 
r : 

M 5 (3M^/3t)^ = constant, independent of r or . (2.14) 

(See preceding paper — Thorne 1976 and §11. e of this paper— for notation used 
here and below. ) Also assuming the star is stable, the stellar structure is 
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(2.15) 


stationary on timescales << 10 


7 


years : 


[-7^7 (any stellar-interior variab le) ] £ ixed r = 0 


This stationarity , together with the identity 

r\ 1 M r\ 

%t^r = %tV + 


■3t 


‘DM 


(2.16) 


and equation (2.14), implies a simple relationship between the time deriv- 
atives and the radial derivatives which appear in the stellar-structure 
equations (T,ll): 


(— -) = - M(— ^— ) 

^3t^M k 3M \ 

r r z 


(2.17) 


For example, if we let both sides of equation (2.17) act on the radius 
function r , and if we combine with equation (T,lla), we obtain the rela- 
tion 


M = 47rr p v . $ 

in 


(2.18) 


where v , the locally measured velocity of inflow of rest mass, is 
m 

defined by 


v in = -(V/R)(3r/Dt) M 


(2.19) 


(cf. eq . [T , 7 ] ) . 

The above considerations are patterned after Paczynski's (1970b) 
analysis of Newtonian stars with mass inflow through stationary shell 
sources. When equation (2.17) is inserted into the relativistic equations 
of stellar structure (T,ll), it produces the relativistic analogue of 
Paczynski's stationary-shell-source equations: 
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(2.20a) 


dr/dM^ = (4TTr 2 p •y ) 1 , 

dM /dM = Zhf , 
tr r 

d$/dM = [GM r /(47rr 4 p) ] Snr 

d(L ft 2 )/dM =« 2 (r — £ ) + RM[dII/dM - (P/y 2 )dp/dM ] , 

r r nuc oV r ' r 

d(L nV & 2 )/dM = a 2 e 
r r nV 


o\J 2 . 2 

d(L r sn/dM r = B e ov 


dX. /dM 
x r 


fffia. /M 
I x 

1 0 


9 


if 


V ^ < V , 

rad ad 


if 


V a > V . 
rad ad 



(2.20b) 

(2.20c) 

(2. 20d) 

(2 . 20e) 

(2. 20f ) 

(2 . 20g) 

( 2 . 20h) 

( 2 . 201 ) 


Here we have replaced all partial derivatives (3/3M ) by ordinary deriva- 
tives d/dM because all time derivatives have disappeared from our equations. 
In equation (2.20g) we have imposed the physical constraint that the abund- 
ances not change radially in the convective region. 

At the outer edge of the middle region the relativistic correction fac- 
tors S, W, ft, If all differ from unity by <10 (cf . eq . [2 . 1] ) ; 

temperatures are so low that no nuclear burning has occurred; and consequently 
the above equations of structure for the middle region reduce to the stand- 
ard Newtonian equations of structure with constant luminosity, which we use 
in our outer-region analysis. This fact guarantees that we obtain a reason- 
able match between middle region and outer region by simply enforcing con- 
tinuity of the fundamental variables r, M , M , $ j L r J > L °’ > X j_ > 

T . and c at radius r . But in doing so we must be careful with M 
3 H o-m r 
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and M : In Newtonian theory M is both rest-mass and active 

tr J r 

gravitational mass. In general relativity is rest mass, 

while M is active gravitational mass; and the additive normalization 
of M is crucial, while that of is unimportant (all details of the 

model except are unchanged by the renormalization constant) . 

These facts dictate that 

M of outer, \ / M fcr of middle, 

Newtonian region/ \Relativis tic region 

one need not enforce any matching condition on the middle-region . 

(2.21b) 

In our analysis of the middle region we use specific analytic expres- 
sions for all the auxiliary variables (pressure P, opacity K , relativistic 
correction factors £, <# , R, 'y , etc.) as functions of our fundamental 
variables : 

The nuclear species which we consider are (hydrogen) with abundance 

/ 12 
X = Xjj, ne (helium with abundance Y = X^, c (carbon) with abundance 

C = X , and "metals" with abundance Z H 1 - X - Y . The binding energies 

per baryon relative to hydrogen are 

helium: (1 - m^Mm^c 2 = 0.007118 c 2 (2.22a) 

carbon: (1 - n^/lZm^c 2 » 0.007118 (l + l/11.0)c 2 , (2.22b) 

products of carbon burning: 0.007118 (1+ 1/11.0 + 1/13) c 2 ; (2.22c) 

and, consequently, the mean binding energy per baryon is 

1 _x_v 1 _v_Y_r o 

B = 0. 007118]; 1 - X + YT7T + M'f “ "1 c ~ * (2.23) 

In our numerical calculations we have assumed that hydrogen burns by the normal 


) at join point r ; (2.21a) 

/ J - o-m 
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CNO cycle; this is a serious source of error, as will be discussed in 

§V. For the normal CNO cycle 93.6 percent of the energy goes into heat 

CNO 

and 6.4 percent into neutrinos; hence, the energy generation rates £ ^ 

and £ C ^ are related to the rates of change of hydrogen and helium abund- 

n\> 

CNO , CNO , 
ance ot^ and 0^ by 


.CNO , _ 2, CNO 

, = (-0.006662c )a , 

nuc H 


CNO = (_ 0 .000456c 2 )a„ NO 
n\> H 


(2.24) 


CNO 

He 


CNO 


CNO 

The CNO energy generation rate £ nuc is expressed as a function of X, Z, 

p, and T by equations (17.280), (17.282), and (17.283) of Cox and Giuli 

(1968) with X = Z/2 . When helium burns by the 3a process to form car- 
UJN 

bon, neutrino losses are negligible; hence 


3a 


2. 3a 


3a 


£ = (-0.000647c )al , = 0 , a' 

nuc He nV 


,3a 


3a 


■“a* 


(2.25) 


We use equations (17-341) and (17.342) of Cox and Giuli (1968) for the 3a 

,3a 


energy generation rate £ 


nuc 


12 12 

We assume that carbon is burned by C + C 


reactions, and in doing so we ignore neutrino energy generation: 


CC , A 2\ CC CG „ fry ryf \ 

£ = (-0.00055c ) a„ , £ = 0 . (2.26) 

nuc c 

We use the Arnett-Truran (1969) analytic expression for the CC burning rate 
together with the Salpeter-Van Horn (1969) analytic expressions for the 
screening factors. The non-nuclear-burning neutrino energy generation rate 
£ QV (X,Y,p,T) (including pair, photo, bremsstrahlung, and plasma neutrinos 
but excluding URCA) we take from Beaudet, Petrosian and Salpeter (1967). 

The pressure P and specific internal energy IT are split up into 
four contributions: radiation, ions, ionization electrons, and pairs: 


P = P 


rad 


+ P. + P + P 
xon e 


it = n 


pair 


rad 


+ 


it. + n 

xon e 


+ n 


paxr 


(2.27) 


22 



The radiation contribution has the usual form 

P . = \ aT 4 , n = 3(P ,/p) , a = 7.5647 X 10~ 15 erg cm“ 3 K“ 4 , 

rad 3 rad rad 

(2.28) 

Crystalization of the ions is ignored in P and IT ; they are assumed to 
form a perfect gas with mean molecular weight 


y . = (X + Y/4 + Z/16) 

ion 


(2.29a) 


for which 


P. = (k/nu) (p/y . )T , n. 
ion H ion ion 

In the middle region temperatures are so high (T 
fully ionized, and the mean molecular weight per 


3 

2 


(P. 

ion 


/P) • 


(2.29b) 


> 10^K) that the plasma is 
ionization electron is 


y e = 2/(1 + X) 


(2.30) 


Our middle region covers a temperature-density regime in which the ioniza- 
tion electrons range from extreme nondegeneracy to extreme degeneracy (see 
Fig. 2). Over the entire range we describe P e (T,p,P e ) and IT(T,p,y e ) by 
the Eggleton-Faulkner-Flannery (1973) analytic fit to the relevant Fermi- 
Dirac integrals; in that fit we use their "thermodynamically consistent 
coefficients" (their Table 5) . Near the knees of our supergiant models 

electron-positron pairs play a crucial role (see Fig. 2 and the discussion 

6 —3 2/3 

in §V) . Fortunately, the pairs occur only in a regime [(p/4xl0 g cm ) 

2 

<< kT/m c << 1] where their contribution to P and IT can be expressed 
e 

in fairly simple analytic form and can be added linearly onto the contribu- 
tions from other sources. The relevant expressions for p , and H . are 

pair pair 

given in the Appendix [eqs.(A.8) and (A. 11) ] . 
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In the middle region all sources of opacity are negligible except 
electron and positron scattering of photons, and opacity to heat conduc- 


tion by degenerate electrons: 


K = (1/k + 1/k , ) 

v es deg.e 


-1 


(2.31) 


We use the following analytic formula for the scattering opacity 


K - (0.4 cm 2 /g)u” 1 (l + 2n,/n ) G(T) 
es e + e 


(2.32a) 


G(T) = 0.4+ 0.6 exp(-0. 04328 T ) , if 0 < < 20 (2.32b) 

G(T) = -0.13887+ 4.9871 T~ 1/2 - 5. 9479. T" 1 - 2.362 T~ 3/2 , (2.32c) 

if 20 < T. < 125. 
k 

Here = kT/(l keV) = T/ (1.160 X 10 7 K) ; G(T) is the special relativistic 
correction to the electron-scattering opacity; formula (2.32c) for G(T) 
is taken from Sampson (1959); formula (2.32b) is our analytic fit to 
Sampson’s computations; and n + /n e is the number density of positrons 
divided by the number density of ionization electrons as given by equations 
(A. 10) and (A. 8) of the appendix. At the time of our numerical work the 
Flowers-Itoh (1975) degenerate-electron heat conductivities were not avail- 
able, so we used Paczynski’s (private communication) analytic fit to the 
tables of Canuto (1970) for carbon: 

log 10 K deg.e = “°* 05 + 0,533 Pg " 1,057 1 o 8 10 P 6 + 2,17 lo 8 10 T 9 * 

(2.33) 

(This formula gives a good fit for 1.05 + 3 log^T^ < log^pg < 6.15 + 

6 3 

3 log^Tg.) Here is density in units of 10 g/cm , and is 
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g 

temperature in units of 10 K. 

_8 

Because the total amount of mass in the middle region 10 M ) is 
negligible compared to that in the core, the gravitational field in the 
middle region is (very nearly) the Schwarzschild gravitational field of 
the core: 

$ = ~ c 2 £n(l - 2 GM /c 2 r) . (2.34) 

2 tc 

Here M is the total mass ("active gravitational mass") of the core 


M = M at outer edge of core. (2.35) 

tc tr 


In our middle-region computations we used expression (2.34) for $ ; we 
used the corresponding Schwarzschild-metric expressions for the redshift 
and volume correction factors 


, 2 1/2 

(1 - 2GM /c r) 
tc 


and we used expressions (T,6c,d,e) for the relativistic correction factors 

& , e , ar . 

Our computation of the radiative, adiabatic, and convective gradients 
V ^ad* anc * ^ l°ll° wec l the prescription of equations (T,9) and (T,10) 
with mixing length equal to pressure scale height = Hp . However, in 
our solution of the mixing equations (T,10) we stupidly used Newtonian 
rather than relativistic expressions for g and H : 


g - GM /r = g .hr = (0 . 84 to 1.0) X g 

tr “correct “correct 


(2.37) 


H = P/pg = (H J x W= (1.00 to 1.43) XH 

p p correct p correct . 

These errors have the same effect on the star's structure as using the cor- 


rect g and H , but making the ratio £ /H increase from its chosen value up 

p t p 
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to 1.25 its chosen value as one moves inward toward the knee of the star — 
i.e., they cause the convection to be a little more efficient than it should 
have been near the knee. Because the mixing-length theory is so unreliable, 
and because the convection is fairly efficient near the knee, we have not 
recomputed our models with these errors corrected. 


e) Global Structure of the Computation 

To construct a 3tellar model one can proceed as follows: (i) Specify 
the following parameters: 


(X , ,Y ,C . ) 
ph ph ph 




M 


t,m-i 


R . 
m-x 


= (photospheric abundances of hydrogen, helium, and carbon) , 

5 (total mass of star) = ("active gravitational mass"), 

8 

= (total mass of inner region) = M tc [l+ an error of 0(10 )]. 

-2 

= (radius of inner region) = R [1+ an error of 0(10 )]. 

(2.38) 


For given M . the value of R .is taken from the theory of bare 
t,m-x m-x 

(envelope-free) neutron stars. (ii) Pick trial values of the quantities 
required for starting inward integrations: 

L = (total photon luminosity of star) , 

Tph = (photospheric temperature) j 

and also pick a trial value of 

M = (rate of inflow of rest mass), 

which plays an important role in the middle region but not the outer region, 
(iii) Integrate the equations of stellar structure inward from the photo- 
sphere to the middle-inner match point r . ; and iterate the three trial 


(2.39a) 

(2.39b) 

(2.39c) 
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parameters L, T^, and M until the three matching conditions 


8 , 3 


M — M . , r = R . , L =0 at p = 3 X 10 g/cm 
tr t,m-i m-i’ r 


(2.40) 


(eqs.2.13) are satisfied. 

In practice the L-T^-M parameter search is not difficult: 

L and T , can be determined with rather good accuracy by Newtonian, 
ph 

outer-region integrations only (see §IV.d) — and these can be performed 
once and for all, with ease, to give a family of outer-region models for 
subsequent join onto middle-region models. Moreover, in the case of giant 
stars, where negligible nuclear burning occurs in the convective region, 
and where— it turns out — non-nuclear-burning neutrino losses are negligible, 
one can express M as an analytic function of M , R , and L . But in 
supergiants hydrogen burning in the convective envelope prevents one from 
finding an analytic expression for M . 

The principal key to the giant-star expression for M is the follow- 
ing equation of energy conservation, which is valid everywhere in our 
stellar models except in convective nuclear burning regions: 

+ L^ V = ft ^ [Mc^ yft + constant] 

= M(H + P/p - B + $) + constant in Newtonian limit 

(2.41) 

[This equation can be derived as follows: (i) add eqs . (2.20d,e,f); 

(ii) use eqs. (T,20) and (2.17) to eliminate £ nuc + £ nv (this step re- 
quires that the nuclear-burning region be non- convective); (iii) write 
-(P/p2) 3p/3M r as 3(P/p)/3M - p ^3P./3M , and use eqs. (2.20i,c) and 
(T,6a) to express 3P'/3M r in terms of 3 ft /3M ; (iv) use definition (T,6e) 
of U to bring the equation into perfect differential form; (v) integrate 
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it.] Another key to the expression for M is a conservation law for 
the nuclear-burning-induced neutrino losses L^ V , again valid everywhere 
except in convective nuclear burning regions: 

L nV = ft 2 [ 0.000456 Me 2 X ft u + constant] . (2.42a) 

r H 

Here ft„ is the value of ft at the center of the hydrogen burning 
n 

shell, which is so near the core that 

ft u = ft = (1-2 GM / c 2 R ) 1/2 (2.42b) 

He tc c 

is a good approximation. [This equation can be derived as follows: 

(i) in equation (2.20e) replace £ by expressions (2.24), (2.25), 

CNO 

and (2.26), and then replace a by expression (2.20g); (ii) invoke 

ri 

the fact that X changes only in the hydrogen-burning shell, which i3 so 
thin that it has ft essentially constant throughout; using this fact 
write the equation in perfect differential form; (iii) integrate it.] In 
giant stars it turns out that the non-nuclear-burning neutrino losses are 
totally negligible throughout the outer and middle regions 
(L°^ - ft 2 L°^ << L ) , and no significant nuclear burning occurs in convec- 
tive regions. Thus, equations (2.41) and 2.42) can be combined to obtain 
the following relation, valid throughout the outer and middle regions: 


L = ft _2 L + ft 2 Mc 2 [£4ft - 1+ B Jc Z + 0.000456 ft„(X , -X) ] . (2.43) 

r ph H ph 


Here the constant has been evaluated at the photosphere, where L^_ = L , 
ft can be approximated as unity, U = 1 - with the photo- 

spheric value of the nuclear binding energy, and is the photospheric 

value of the hydrogen abundance. To obtain the desired expression for M, 
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If 


we need only evaluate expression (2.43) at the inner edge of the carbon-burning 

o 1/2 

shell, where L = 0,X = Y = C = 0, ft = (1-2 GM /c R ) , and U 

r t c c 

2 

can be approximated as 1 - B/c with B taken from equation (2.23): 

• 2 f / 2GM tc U/2 

Me = L -( 1 — - (0.991687 + 0.000456 X , ) + 1 - 

L ' c 2 r c > P h 

/ 1 - X - Y u l-X.-Y.-C.v 

- 0.007118 (l- X ph+ t + *-) 


-1 


(2.44) 


III. NUMERICAL MODELS 

Some details of our numerical models for stars with neutron cores are 
shown in Figure 2 and Tables 1-4. The physical features of these models 
will be discussed in §§IV and V. 


IV. DETAILS OF THE STELLAR STRUCTURE: GIANT MODELS 

Table 1 and Figure 2 display the internal structure of a typical 

giant model — one with a total mass of 5 M and core mass and radius of 1 M 

0 o 

and 10 km. Tables 3 and 4 show some details of other giant models. In 
this section we shall point out and analyze analytically some important fea- 
tures of these models. 

a) Overall Structure 

In §IIb we explained, analytically, the hydrostatic and thermal de- 
coupling of the core (inner region) from the rest of the star. We shall 
now elucidate the reasons for the gross features of the rest of the star 
(extremely thin halo surrounded by very deeply convective envelope) . 

Consider the forces which act on the plasma (gas) in the nondegenerate 
region r > R c . The pull of gravity is counteracted by the plasma's own 
pressure-buoyancy force and by the force of outflowing radiation: 
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SKIT 

g 


(4.1) 


gravitational force\ 
per unit volume 


gm p 

r 

2 

r 


T rad 

- dP KpL 
1 _g. r 


V dr 


. 2 
47Tr c 


Here 

only 


V 

g 


is the relativistic 


enthalpy correction factor for the gas 


U = 1 + (n - B + P /p)/c 2 (4.2) 

g g g 

(cf. eq.[T,6e]), L is the locally measured luminosity carried by dif- 
fusing radiation 


L rad = 
r 


L 

r 


4ttA 


F 

conv 


9 


(4.3) 


and all other quantities have been defined earlier (cf. §I.e). This force- 
balance equation can be derived either from first principles, or from the 
relativistic equations of stellar structure (2 .20a, h, i) , (2 .28) , (4. 2) , (4. 3) , 
(T,6d,e) , (T,8a,b) , (T,9a) , and (T,10a). By analogy with Newtonian theory, it 
is convenient to rearrange the force-balance equation (4.1) as follows: 


where 



'If dr 


GM p 
r K 

2 

r 


&K ir 

g 




rad 

i-/ 

r 


crit 

r 


pri t" — 1 

L = 4ttG cM k <&u It 

X X ^g 


(4.4a) 


( 4 . 4b ) 


is the ''critical luminosity" above which the force of outflowing radiation 
on the plasma exceeds the force of gravity. Equation (4.4a) shows that 
the scale height for the gas pressure is 
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H 


-‘Zrdr 
’g d in P 


,P wGM .-1 
'pc' ' 'c -r 


g 


T rad 
1 (l 

\ cnt / 


-1 


" 1 <4 - 5) 


, rad , 
L \ -I 


(Here use is made of the plasma equation of state P^ = (p/yrn^) kT. ) 

As one moves outward through the star, this equation for Hp first 

6 

becomes valid where electron degeneracy turns off (at p ^ 10 g/cm , 

r = R^) . At that point all quantities on the right-hand side of the 

equation are of order unity, so H p /r 'v 6 X 10 ^ . Thus, the plasma just 

g 

above the core's edge has the extremely small scale height of a hot 

g 

(T = 10 K) neutron-star atmosphere: Hp ^ 1 meter. Physically this scale 

g 

height is governed by the inability of the mean particle kinetic energies, 

-4 2 

kT 'V 10 m^c , to compete with the extremely strong pull of gravity, 

GM /c 2 r ^ 0.1. 
c 

The "halo" of our models is the region just above the core where 
Hp ^ 1 meter. As one moves outward through the halo a distance 15 meters, 
the density drops to ^ 10^ x e ^ ^ 1 g/cm^. This rapid density drop cannot 
continue for many more meters if the star is to support a massive envelope 

around itself. Something must happen soon to increase Hp /r from 

_4 g 

^ 6 x 10 to ^ 1 . Equation (4.5) shows two ways to increase Hp /r : 

-3 g 

(i) by a decrease of the mean molecular weight to a- 10 due to a profuse 

■v “p\r\ f'T* *1 

turn-on of electron-positron pairs; (ii) by an increase of /L^ to 

near unity so that the force of outflowing radiation on the plasma strongly 

counteracts the inward force of gravity. In all of our models the radiation 

force (case ii) is responsible for the increase in Hp /r . It is con- 

g 

ceivable — but seems unlikely to us — that one could build models of type (i) , 
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where H /r increases due to profuse pairs . 

g 

In our giant models , as one moves outward through the halo (where 

energy transport is all radiative) , gravitational energy release drives 

l-cid 3 

= L up higher and higher. Ultimately, at p 1 g/cm , the lumin- 
osity goes supercritical and /r becomes 'v 1 . Very shortly 

g 

before this point the force of outflowing radiation on the plasma becomes 

1 

so great that it begins to drive convection. The onset of convection 

^The Newtonian proof (Joss, Salpeter, and Ostriker 1973), that convection 
sets in before becomes supercritical, is easily generalized to rela- 
tivity theory. 


marks the end of the halo and the beginning of the convective envelope. 

Throughout the strong-gravity region of the convective envelope, the 

plasma is protected against the pull of gravity by the force of outflow- 

ra .d crit -3 » 

ing radiation (1 - /L ^ 10 ). Because the radiative luminosity 

is so extremely close to critical, the star is forced to remain convective 

throughout this region. Ultimately, with increasing radius, gravity 

weakens enough that there might be some hope of the plasma supporting 

itself without the help of radiation forces. However, the outflowing 

luminosity cannot now be shut off. It is pouring outward with a rate L 

exit 

designed to counterbalance gravity at small radii, L^_ - L_^ (strong- 

gravity region) ; and with ever-increasing r and ever-decreasing T the 

exit 

opacity is rising higher and higher, driving L r lower and lower. Thus, 
the star remains supercritical (and therefore convective) all the way from 
its knee out to the photosphere. 
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One knows from the theory of stellar envelopes that because our 
stars have very deep convection they must be near the Hayashi track of 
the H-R diagram where photospheric temperatures are low 

T , < 300 OK . (4.6a) 

ph ~ 

The above argument shows, moreover, that the luminosities of our stars 
must be 

L 'u L cr ^ t: (s trong-gravity region) ^ 4ttGcM /k 'u 4- x 10 . (4.6b) 

These numbers agree with the detailed models of Tables 1-4. 

We suspect, but are not certain, that it is impossible to construct 
equilibrium models of stars with neutron cores and massive envelopes 
that lie elsewhere in the H-R diagram. The extreme force of gravity near 
the core probably always enforces deep connection and very high L — and 
thereby red-supergiant surface features. 

Using the above information about the stellar structure, we can un- 
derstand semiquantitatively the flow of energy inside the star: Mass 

flows from the static envelope, through the inflowing envelope, into the 
halo, and thence into the core. In the inflowing envelope, because of 
inefficiency of convection, the temperature gradient is slightly super- 
adiabatic, so the inflowing matter gets heated not only by adiabatic com- 
pression due to gravity, but also by the absorption of some of the up- 
flowing luminosity . Mathematically, in the equation of energy 

generation (2.20d) dll/dM^ - (P/p )(dp/dM ) is negative due to super- 

2 

adiabaticity , so LB (the redshifted luminosity) increases inward. 
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Equivalently, in the equation of energy conservation (2.43) superadiab- 

. . 2 
atxcxty means that increases inward, so L fi also increases inward. 

r 

By the time it reaches the knee, the inflowing matter contains an 
enormous amount of internal energy, almost all of it tied up in radiation: 

(wa) K > (&») h . — -» 

n . ii+ p/p . 

____ * ■ .> a _ B / c 2)[(l- 2GM /c R ) 1/Z - 1] 'U 0.2 . (4.7) 

c c Ph L *' 

[Here we have used expressions (T,6e) and (2.36) for %C and R, , together 
with the fact that because the knee is so close to the core boundary, the 
redshift factor B is very nearly the same at the knee as at the core 
boundary. ] At the knee the temperature gradient goes very sub adiabatic 
(in fact, nearly isothermal) , so the contracting matter begins to release 
its huge store of thermal energy, converting it into outflowing radiation. 
Because its temperature is now remaining constant, its specific internal 
energy II - aT^/p falls off as 1/p . After the density has increased 
by only a factor 10 above P ^ nee 5 90 percent of the stored energy has been 
converted into luminosity . After several more decades of density 
increase nuclear burning begins to occur, producing further luminosity 
(but much less than was produced by gravity and released just below the 

knee). By the time the flowing matter gets inside the core, essentially 

r • 

all the star's luminosity has been accounted for; has dropped nearly 
to zero. Overall, the contribution of gravitational contraction to the 
total luminosity of the star is 
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(4.8a) 


L grav= - B ph /cZ>[1 - (1 ‘ 2GM tc /R c c2 > 1/2 I 

(cf . eq. [2. 43] and associated discussion); and the contribution of nuclear 
burning is 


. o ') 1 In 

L = Me (1 - 2GM /R c ) ' 0. 007118[0. 936X . + (X Y .)/ 11.0 

nuc tc c L ph ph ph 


+ <V V V /13] • <4 - 8b) 


The ratio L /L = L / (L + L ) is shown for various models in 

nuc nuc nuc grav 

Tables 3 and 4. 

For a detailed example of these features of energy flow, see the 

2 

columns labeled r-r, p , and 1-B L /L in Table 1. 

K r 

Non-nuclear-burning neutrino losses are totally negligible (<< 1 L ) 

8 3 

in the outer and middle regions (p < 3 x 10 g/cin ) of all our models; cf. 
Tables 1 and 2. We have not made a thorough search for models with high 
neutrino losses; but we suspect that high losses are incompatible with 
stellar equilibrium as well as stability. 

b) Structure of the Halo and Sharpness of the Knee 

The halos of our giant models are remarkably isothermal, and the 
transition through the knee into a superadiabatic temperature gradient 
is remarkably sharp (see Fig. 2). These features can be understood as 
follows : 

To avoid issues cf radially changing chemical composition, consider 
that region of the halo which lies outside the hydrogen-burning shell 
(p < 3 x 10 g/cm for the 5 model of Fig. 2), Here the pressure and 
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internal energy due to gas and radiation are 


p , m. n - • 

g Um ’ g 2 p > 
H 

v-l 


14 

P = ^ aT ; n , 
rad 3 rad 




3P 


rad 


= const. (4.9a) 


Because the halo is so thin in radius and contains so little mass, through- 
out it we can set r = R , M = 11 ; and thence 

c 9 tr tc 


r 1 = R = R c = (1 - 2GM tc /c 2 R c ) 1/2 


(4.9b) 


2 3 

(cf. eq. [2. 36] ) . Also, because P/c < p « M tc /^' rTR c throughout the 

halo, and because nuclear binding energies and particle kinetic energies 

2 

are small compared to n^c , we can approximate 


= M /M , % = 1 

tc r g 


(4.9c) 


(cf. eqs.[T,6c] and [4.2]). Finally, because all luminosity is carried 
radiatively in the halo, we can set . 

By using the above relations we can rewrite the force-balance equa- 
tion (4.4a) for the plasma in the halo as 


dP 

dr 


R 




(4.10a) 


where .is the acceleration of gravity at the edge of the core 


g = (GM /R 2 ) ft _1 . (4.10b) 

- c tc c c 

The analogous equation of force balance for the radiation is obtained by 
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setting P ra( j = P - P , by taking the difference of equations (2.20i,a) and 
(4.10a), and by invoking the relations (4.9b, c), (4.10b), and ft" - 1 — 

4 p rac j/p c [cf* eqs. (T,6e) and (4.9a, c)]: 


dP 


rad 


4g P , 
& c rad 


g. 


dr 


£ 2 
c c 


£ T crit 
c L 

r 


(4.11) 


The first term on the right-hand side is a gravitational redshift term; 

it can be neglected because of the thinness of the halo (r- R << R ) : 
° K c c 


dP ,/dr = -(g /£ ) p (L /L Crit ) 
rad °c c r r 


(4.11’) 


By taking the ratio of the force-balance equations (4.10a) and (4. 11’) and 
combining with the equation of state (4.9a), we obtain 


d &n T / d £,n p = (4y /y - 1) 

§ 


-1 


(4.12a) 


where 

y T = L crit /L - 1 = 3 t /( 1-3 t ) , y = p /p , = 3 /(i-$ ). 

L r r L L g g rad g g 

(4.12b) 

We shall see below that y^ » y^ throughout the halo, except very near 
the knee and near the nuclear burning shells; thus, the halo must be nearly 
isothermal (d Jin T / d £n p << 1) . 

The opacity in the halo is due, almost entirely, to electron scatter- 
ing and thus depends on temperature but not density (eq. 2.32)— and is 

crit 

essentially constant throughout the halo. Thus, (eq. 4.4b) is also 

essentially constant with value 


4ttGc M 3.2x10 L o 

crxt tc e tc 2 


£ K 
c es 


£ G(T ) M 1 + X * 
c K ® ph 


(4.13) 
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) 


where G(Tg) I s t^ e Klein-Nishina correction function for the electron 
scattering opacity, evaluated at the temperature of the knee T R . The 
knee occurs where the temperature gradient becomes adiabatic; thus 


d in T / d In p = (9 £nl/3 In p)^ = 1/3 at knee . (4.14a) 

(Here we have used the fact that P , » P near the knee, so that the 

rad g • 

3 

adiabats are T /p = const.) Equations (4.14a) and (4.12a) show that 


Y T = Y at knee ; y T > Y in halo. 

L 'g 'L 'g 


(4.14b) 


The value of Y^ in the halo is governed by , which is determined by 

the equation of energy conservation: 


L = L nUC + ft -1 M(n ,+ P ,/p) = L nUC + const/p . 
r r c rad rad r 


(4.15) 


[See eqs.(2.43), (4.8b), and (T,6e) specialized to the case r - << R^; 

B - B , , X = X , , and II + P /p radially constant because of isotherm- 
ph ph gas gas 

ality. Here L^ UC is the total contribution of nuclear burning to in 

the halo 

L nUC = ft _2 L = fc _1 Me 2 • 0.007188[0.936 X , + (X + Y , )/ll.O 
r c nuc c ph ph ph 


+ < W + V /13] ; 


(4.16) 


cf. eq.(4.8b).] By combining equations (4.15) and (4.12b), using (4.14b) to 

exit 4 

evaluate the constant in (4.15), using theconstancy of , and ignoring a 

tor Y ir where it is unimportant, we obtain 
(P /P K )( ^K + » 


= 


1+8 (p/p v - 1) 

nuc ^ K K 


-1,3 


nuc crxt 
r 


— ^ 0.03 . 


(4.17) 


Here p T , and y _ _ are the values of p and Y at the knee, and 
K gK K g ’ 


fac 
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8 <0.03 because 

p nuc ~ 


c nuc 


L nuc = R" 2 L 'V 0.030 ^" 2 L ^ 0.030 L (knee) - 0.030 L 

c 


crit 

r 


Because y = y (p/p ) in the isothermal region, we have 
g gK K 


Y L , (P/P K 1)(1 ' ^nuc" 0 nuc Y gK p/( V 
— = 1 + 

Y g [1+3. 


nuc (p/p K“ 1)](Y gK p/p K* 


(4.18) 


The isothermal region is that region in which y /y >:> Y (cf. eq. [4.12a]). 

o 

Equation (4.18) shows that it extends over the range 


5 x 10 ~ * Y gK « (P/P K - 1) « ^^nucV ^ 7 X 10 


+4 


(4.19) 


The p-T curve for the 5 M star in Figure 2 demonstrates this: At the 

© 

left end of the halo the termination of isothermality is so sharp 

( Ap/p a. y ^ 5xl0 _Zt ) that the slope of the p-T curve looks discontinu- 
gK. 

4 3 3 

ous. Toward the right isothermality ends at p 'v 2 x 10 p^^lO g/cm . The 
above analysis diagnoses correctly small departures from isothermality; but 
as the departures become significant (d £n T/3 &n p ^ 0.1), the analysis 

produces serious errors. 

The density-radius relation in the isothermal region can be derived 
by combining the plasma equation of state (4.9a) and expressions (4 . 12b) , 
(4.17) for L r /L^ rit: with the plasma force-balance equation (4.10a), and 


then integrating. The result is 

p - p., = constant x exp 1 


vvv 


(4.20a) 


where is the value of the gas-pressure scale height (4.5) a few meters 


below the knee where L = L^ UC 
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(4.20b) 


For the 5 M star of Table 1 and Figure 1 this formula gives ft H, =1,30 
0 C “ 

meters. The density profile (4.20a) agrees rather well with the numerical 
model of Table 1 inside its realm of validity (eq. [4.19]). For example, it 
describes within a few percent accuracy the increase in density scale height 

from 11^ = = 1.55 meters deep in the halo to =H^(1 - p^/ p) ^ =50 meters 

—9 

at 1 - p^/p = 3 x 10~ , However, very near the knee (at 1 - P R / P ^ Y gK = 

5 x i0 _/t ) , it breaks down because of the breakdown in isothermality 
(d &n T / d ftn p no longer « 1). 

Unfortunately, in the neighborhood of the knee there is a serious 
omission in the physics which we have put into our analysis: We have ig- 

nored the possibility of "convective overshoot" in which turbulent cells 
plow through the knee and into the upper layers of the halo before being 
stopped by pressure buoyancy forces. 

We can estimate the effects of convective overshoot in our 5 M model 
(Table 1) as follows: Just above the knee the mixing length (assumed equal 
to pressure scale height) is - R^/4 = 2.5 km (cf. eq.[4.26] below and 

recall that P « T^) . Table 1 shows that convective cells within this 

7 8 

distance of the knee have typical velocities of v 10 cm/sec. Suppose 
that a small cell moving downward with this speed hits the knee, and that 
when it hits it has the same density and temperature p R and T^ , as its 
surroundings. Because the cell's velocity v^_ is far less than its sound 
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9 . 85 

speed (v - 10 cm/sec), it maintains pressure equilibrium with its 
s 

surroundings as it penetrates the halo. Pressure equilibrium means tem- 
perature equilibrium since P , » P , which means constancy of tern- 

rsu gas 

perature since the halo is isothermal. Assuming negligible heat exchange 

3 

between the cell and its surroundings (T /p constant in cell) , we con- 
clude that the cell maintains constant density, P^-j^ = Pjr > as it pene- 
trates the halo. Consequently, its deceleration rate as it moves 
through the halo is given by 

p^ 1 dv / dt = " 8 c ( p “ P K > 

where g^ is the (constant) acceleration of gravity throughout the halo 

_2 

(eq. [4.10b]). Since v = - dr/dt , and since the density profile has 
the form (4.20a), we can rewrite this deceleration equation in the form 


dv /d(p- P K ) = - 2g c H h /p K 


(4.21) 


Integrating this equation and imposing the boundary condition v = v at 
p = , we obtain for the density P 0vers | l00t at which the cell halts 

its plunge and begins to rise 

2 2 2 


v t u ”h c 


Vers h oot _ . -1(1- 6 ) . (4.22) 

P K 2g cHh 2 l n»c' c 2 M k 


7 8 

For our 5 M model, with .3 = 0.028, v = 10 cm/sec, y = 0.62, and 

a ’ nuc t 

8 2.3 

T K = 10 * 'K this gives P overs ^ oot /P^ ~ 1*08 . Cells moving three times 
as fast will penetrate 10 times farther, i.e., to p , /p Tr - 2 . 

The above estimates suggest that convective overshoot is of some, 
but not great importance. However, the following factors make this con- 
clusion somewhat uncertain: (i) We evaluated the convective overshoot 
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assuming a small convective cell, but the size of a typical cell just above 
the knee is probably Ml^/10 ^ 1 km, which is far greater than the scale 
height (a few tens of meters) of the region into which the cell penetrates, 
(ii) The region of overshoot is the region of greatest gravitational energy 
release: the energy release between p and p is 

AL r = (L^ rit .- L™ C )(1 - p K /p) ; (4.23) 

cf. eq.(4.15). A serious modification of the temperature distribution in 

this region due to overshoot will seriously affect the details of gravita- 
release, 

tional energy/ and will thereby affect the average density profile and the 
pressure-buoyancy force on the convective cell, and might thus seriously 
affect our above estimates. 

Obviously, a detailed study of overshoot is needed. 

In this discussion of the halo, turn attention now to the nuclear 
burning shells. Because of the extremely small scale height in the halo and 
in the outer layers of the core, the nuclear burning shells are very thin: 
typically (physical thickness) = B ^ Ar 'v. 2 meters for hydrogen shell, 

4 meters for helium shell, and 20 meters for carbon shell (see Table 1), 

The time required for matter to contract through these shells is 
B ^Ar/v. o, 10 sec for the hydrogen shell, 'V 30 minutes for the helium shell, 
and % 10 days for the carbon shell. Note that the electrons are nondegener- 
ate in the hydrogen shell, slightly degenerate in the helium shell, and 
fully degenerate in the carbon shell. However, these conclusions, particu- 
larly concerning the hydrogen shell, are somewhat uncertain because of 
inadequacy of our nuclear burning rates (cf. §V) . On the other hand, 

L /L - 0.030 is so small that errors in our treatment of nuclear burning 
are probably unimportant for the overall structure of the star. 
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c) Structure of the Inflowing Envelope 


In the inner regions of the inflowing envelopes of our giant stars 
7 

(10 K < T < T ) convection is fairly efficient, so the temperature gradient 
K. 

is not far from adiabatic; cf. Table 1, where adiabaticity would mean con- 
stancy of 3 , and Figure 2 where adiabaticity would mean a T-p curve 

o 

parallel to the 3 - constant lines, (One must not diagnose adiabaticity 

S 

from V - V a( j in regions where 3g « 1 ♦) 

By approximating the temperature gradient as adiabatic, we can derive 

simple expressions for the structure of the inflowing envelope. Adiabaticity 

of the flow implies (by virtue of the relativistic Bernoulli equation, or 

2 

by eq.[2.43] with R constant) that is independent of radius. 

Because 3 << 1 and because nuclear binding energies can be ignored, equa- 

g 

/ 2 * 

tion (T,6e) for U reduces to 1 + (4aTj/3pc . By combining this expres- 
sion for U with the relation 


6 g = P g /P rad * (3k/u V )(p/T 3 ) 

and with expression (2.36) for R , and by setting &R = (&R) K , we 
obtain 


g yny 


kT . 4 kT K 

„2 ^ + 3 2 


g y^c 


\/ l - 2GK tc /c 2 R \ p / 2 
7V - 2GM /c 2 r } 


- 1 


(4.24) 


In order that the temperature T not go negative and not go isothermal at 


r » r , the knee temperature must satisfy 
K 


kT. 


K 


1 


(1- 2GM /c 2 R ) 1/2 
tc c 


(4.25) 


gK ym H° 

Deviations from this relation are a measure of the deviation from 


adiabaticity. For the 5 model of Table 1 this relation predicts 
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log T = 8.217 compared to an actual value of log T = 8.249. Using 
relation (4.25) we can rewrite equation (4.24) for the temperature profile 


as 


4_ kT 

8 

g 


GM. 


(1- 2GM tc /c 2 r) 1 ^ 2 


- 1 ^ 


tc 


2 

c r 


(4.26) 


where is the Newtonian limit. 

3 

Note that T “ 1/r implies p 1/r , which means that and 

M increase only logarithmically with radius 


M tr - M tK “ 47r P K R c * n(r/, k> 


(4.27) 


This accounts for the very small amount of mass contained in the inflowing 
envelope (third column of Table 1) . 

d) Structure of the Outer Region 

The outer regions of our models (r > r Q static envelope, photo- 
sphere, atmosphere) are very similar to the outer regions of red supergiants 
with white-dwarf cores or nondegenerate cores. Therefore, we shall not 

comment on their detailed structures or on their sensitivity to the choice 
of mixing length (Table 3) , 

However, it is very important to notice that the luminosities and 

photospheric temperatures, L and T are exceedingly insensitive to the 

details of the core, halo, and inflowing envelope. L and T , are 

ph 

fixed almost completely by the total mass M , the core mass M . and 

t tc 

the envelope composition X ph , Y^ h , C ph . Compare, for example, the fol- 
lowing three models with the same U , M. , X ,, Y 1 , C , , and ratio of 

t tc ph’ ph’ ph’ 

mixing length to pressure scale height: the third model in Table 3 

and the third model in Table 4 (relativistic models with different core 
radii) , and the eleventh model in Table 3 (which is Newtonian) . Despite 
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the difference in their inflowing envelopes, halos, and cores, their 
luminosities and photospheric temperatures agree almost exactly. 

Figure 3 explains this remarkable fact. Figure 3 is an H-R 
diagram for static stellar envelopes near the Hayashi 

track of a 5 M star. All the curves in Figure 3 were constructed using 
e 

* 

Paczynski’s computer program GOB for static stellar envelopes (§11. c) , with 
no attempt to join the envelopes onto any kind of core. Notice the ex- 
tremely narrow range of photospheric temperatures on the horizontal axis . 

The envelopes of Figure 3 can be joined onto a variety of types of 
cores. In the case of a white-dwarf core with hydrogen-burning shell 
source, the base of the static envelope, r Q _ m of eq.(2.1), is near or in- 
side the shell source: thus log T ^ 7 to 8 and stars with white-dwarf 

o-m 

cores typically lie between the solid curves 8 and 7 of Figure 3. 

In the case of a neutron core, the temperature falls off roughly as 

3 9 

1/r between r Tr and r ^10 r Tr ; and because T_ < 10 K, we must have 
K o-m K K. 

6 

T <10 K» In fact, all of our detailed giant models (Tables 3 and 4) 
o-m ~ 

have 5.9 < log T q m < 6.4 . In the envelope H-R diagram (Fig. 3) such 

models lie along an extremely narrow strip, A log T ^ - 0.001 ; and for 

given core mass M , the luminosity within this strip varies by only 
tc 

A log L 'v- 0.02. Thus, to within A log T ^ - 0.001 and A log L 'v 0.02, 
the envelope is oblivious of the details of the core and halo. 

This behavior is due to the well-known fact that as one moves right- 
ward in the H-R diagram, approaching the Hayashi forbidden region, the 
characteristics of the base of the envelope change extremely rapidly. 

The above discussion shows that, for given L , the photospheric tem- 
perature is not even sensitive to the difference between a white-dwarf core 
and a neutron core. The star with neutron core will be redder by only 


A log T ph ^ 0.01. 
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V. DETAILS OF THE STELLAR STRUCTURE: SUPERGIANTS AND MASS GAP 

Consider a sequence of models with fixed core properties (M tc »R ) an ^ 
successively higher total mass M — e.g., the sequence in Table 3. The 
low-mass models have "giant" structures of the type discussed in §IV. The 
high-mass models have "supergiant" structures (convective envelope dips 
into hydrogen-burning shell, and most of energy generation is by hydrogen 
burning rather than gravitational contraction). Between the giant and 
supergiant models there is a "mass gap" in which our computations have 
failed to produce any equilibrium configurations. 

This peculiar situation can be understood as follows (see Fig. 4). 

C1TX t 

The critical luminosity L^ in the inflowing envelope has the form 

cnt = /f7TGc M /k (5.1) 

r tc es 

(eqs.[4.4b], [4.9b,c]), where K eg is the electron-scattering opacity 

K = (0.4 cm 2 /g)[(l + X )/2](l+2n /n ) G(T) (5.2) 

es ph t* e 

cir it 

(eq. 2.32). The product IftL^ is plotted, as a function of temperature 

T for 0.70, M tc = 1 , and p ^ (10 g/cm^) (T/10^K) 2 , in Figure 

4. (The dependence on p , which is exceedingly weak and can be ignored, 

enters through the ratio n + / n e ° f pairs to ionization electrons; see the 

Appendix.) At low temperatures (T < 10 2 K) , R,L^ rit is constant; but at 

T > 10 2 K the Klein-Nishina corrections G(T) begin to reduce the electron- 

cri t 

scattering opacity, and thereby increase fjL^ . At log T =; 8.70 , when 
exit 

RL r has increased by a factor 2.0, electron-positron pairs turn on, in- 
creasing the number of photon scatterers, thereby increasing k , and 

ss 

exit 

thence decreasing RL^ . The turn-on of pairs with increasing T is 
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pr •{ t- 

so sharp above log T = 8.70 that ftL r 
4). 


plummets dramatically (see Fig. 


In the envelopes of our models the local luminosity is everywhere 

supercritical (see §IVa). Moving inward through the envelope, one reaches the 
knee (termination of convection) immediately after goes subcritical. Figure 4 

shows two L^CT) curves, one for the interior of a giant model; the other 
for the interior of a supergiant. The difference between the two is ob- 
vious: The giant goes subcritical, with increasing T , before the peak 

crit . 

of is reached. The supergiant has such a high luminosity that it 

passes over the peak; but shortly thereafter hydrogen burning turns on, 

crit 

driving down through the now plummeting curve. The hydrogen 

burning has to generate a very large luminosity (L nuc - L) in order for L r 


to catch up with the rapid plummet of L 


crit 


The sharpness of the pair turn-on at log T = 8.70 (the sharpness of 
crit 

the peak in ) is responsible for the mass gap between our giant and 

supergiant models. For a model in the mass gap one can choose a total lumin- 

cr it 

osity L such that goes subcritical very slightly before the 

peak (giant-type structure); but such a choice always leads to a knee radius 

r^ larger than the desired core radius R^ — and thus to no viable model. 

crit 

If one chooses L slightly larger, so that L skims over the peak 

and somewhat later plummets due to hydrogen burning (supergiant-type struc- 


ture) , one obtains a knee radius r^ smaller than R c — and again no viable 
model. No choice of L can produce the desired knee radius. 



Unfortunately, the above discussion is based on an inadequate 

treatment of hydrogen burning: Our detailed models utilized a CNO-cycle 

burning rate appropriate to the temperatures of normal stars 

(T 'v. [2 to 10] x 10 ^K) , whereas in our supergiants the hydrogen-burning 
8 9 

shell has T - 10 K. The "hot CNO-Ne cycle" burning rates of Audouze, 
Truran, and Zimmerman (1973) would be more appropriate. However, even 
they would be extremely inadequate: Some of the crucial 8 decays involved 
in the hot CNO-Ne cycle have lifetimes of ^1 to 100 seconds, whereas con- 
vection circulates matter into and back out of our hydrogen-burning shell 
in a time At 'V 0.01 second (cf. Table 2). For this reason we expect 
hydrogen burning to proceed in the following very unconventional manner: 
Convection circulates unburned matter into the hydrogen-burning shell, 
where all strong interactions go to completion almost instantaneously 
(At << 0.01 second). The reaction chains then get hung up waiting for 
3 decays to proceed. After 'MD.Ol seconds the 3 -hung-up matter gets swept 
back up to larger radii (lower temperatures) , where it convectively random- 
walks from place to place, while undergoing 0 decay. Sometime later, after 
the 0 decay is partially or fully complete, the matter random-walks its 
way back into the hydrogen-burning shell, where its strong interactions 
proceed once again. 

In a subsequent paper we hope to analyze supergiant hydrogen burning 
from this point of view. We presume that the reaction products will in- 
clude very peculiar relative abundances of various catalyst isotopes, and 
that these may provide an observational handle for stars with neutron 
cores (see §1 . c) , 

It is quite possible that an improved treatment of hydrogen burning 
will change the hydrogen-shell structure of our supergiants substantially, 
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and will destroy the mass gap between giants and supergiants. 

VI. STABILITY OF OUR MODELS 

We have worried about five possible instabilities in our stellar models: 

Dynamical instability of the envelope , caused by the low adiabatic index 
(T 1 < 4/3) in the regions of hydrogen and helium ionization, where much of 
the envelope mass resides. The situation here is similar to that in red 
supergiants with degenerate white-dwarf cores (cf. Paczynski and ZioXkowski 
1968), since the envelopes there and here are nearly identical. In such 
envelopes the thermal and hydrodynamical time scales are comparable, so 
energy transport has a strong influence on the time development of any in- 
stability. We have analyzed the stability of our envelopes ignoring energy 
transport (stability against linearized adiabatic , radial perturbations); 
see last column of Tables 3 and 4. For envelope masses M < 2 M our 
envelopes are adiabatically unstable; for M v > 2 , they are adiabat- 

ically stable. This result suggests (see, e.g., Keeley 1970a, b; 1975) that 
a more correct, nonadiabatic analysis may reveal either pulsational or dis- 
ruptive instabilities for our least massive envelopes; but that our most 
massive envelopes might be stable against all perturbations, except convec- 
tive ones . 

Thermal instability of the, shell sources . Consider a nonconvective 
shell source with average luminosity and temperature L^ and T , and with 
luminosity and temperature drop across itself of AL^ and At . A crude 
analysis (simple generalization of page 857 of Schwarzschild and Harm 1965) 
shows that an average temperature rise of <5T inside the shell produces the 
following rate of increase of the shell's internal energy: 
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6T 

T 


( 6 . 1 ) 


dE 

dt 


internal 


(AL ) 
r 


v - 8 


AL AT 
r 


Here v is the temperature exponent of the nuclear burning rate, e ^T^. 

The nonconvective halos of our models are extremely isothermal — so iso- 
thermal that SCL^/AL^) (T/AT) has values of ^30 to 40 for the giant model 
of Table 1, and > 1000 for the supergiant pf Table 2. This is sufficiently 
large compared to V that our nonconvective haloes are probably stable 
against thermal runaway (positive 6T sets up a heat flow out of the 
shell which exceeds the increased nuclear burning). Even if the nonconvec- 
tive shell sources turn out to be unstable, their very small contribution 
to the star's total luminosity, and their location deep below the envelope , 
and the very short timescale of the instability will probably prevent the 
instability from producing observable effects at the photosphere. 

In our supergiant models the convective hydrogen shell source should 
be protected from thermal runaway by the 8-decay hangup discussed in §V . 

On the other hand, the convective hydrogen burning described in §V might 
proceed in a series of local flashes.' rather than as a smooth energy flow. 
Even if this is the case, the timescale of the flashes will be so short 
(At << 1 second) that their effects presumably will be smoothed out in the 
overlying envelope. 

Instability of the region of gravitational energy release (p < p < 
lOPj,) . We do not now have any insight into the stability of this region. 
Any adequate analysis would have to take account of convective overshoot. 

Runaway neutrino losses, accompanied by an ever-increasing rate of 
envelope contraction . The computations of Zel'dovich, Ivanova, and 
Nadyozhin (1972) suggest that such an instability may occur in models with 
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halo temperatures much higher than ours — if such models can exist. However, 

a 

our low halo temperatures (T < 1 x 10 K) keep the middle-region neutrino losses 
small (« 1 L q ) and presumably will prevent them from running away. Because 
of thermal decoupling (§11. b), neutrino losses in the core cannot produce 
an instability in the overlying halo and envelope. 

Instability of the mass inflow p attern . Bisnovatyi-Kogan (private 
communication) has argued that the inflowing envelope, halo, and outer core 

may be unstable against perturbations which break the radial constancy of 

• 

M. A specific example of such an instability is the possibility (Cameron, 
private communication) that at densities ^3 x 10^ to ^lO^g/cm^ rapid pycno- 
nuclear reactions and electron capture, followed by intensive neutrino- 
antineutrino emission^ might produce a rapid shrinkage of the outer core. We 
doubt that such instabilities exist, but we have no proof. 

A Henyey-type evolutionary calculation would be a powerful tool to 
use in testing for the above instabilities and others. 

VII. EVOLUTION OF OUR MODELS 

We saw in §IV.d that stars with neutron cores can occupy only an ex- 
tremely narrow strip in the H-R diagram, sitting precisely on the edge of 
the Hayashi forbidden region. The boundaries of this strip can be found 
with good accuracy by static-envelope integrations; see §IV.d and Figure 3. 

Take a star with a neutron core and a given total, mass M^_ , and as- 
sume that it does not undergo any violent instabilities during the time 
required for its core mass M ^ to grow significantly. Such a star should 
evolve through a sequence of quasi-equilibrium states of the type discussed 
in this paper. Restrict attention to giant-type stars, for which the 
envelope does not evolve chemically. Then the evolution of the surface 

features, L and T , can be read off a static-envelope H-R diagram such as 
ph ' ~ 

Figure 3, without any reference to the structure of the inflowing envelope 
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must move up the narrow allowed 


or halo. As increases, L and T , 

tc 5 ph 

strip (strip with T q ^ 10^). The evolution will terminate by collapse 

of the core to form a black hole when M^_ c reaches the Oppenheimer- 

Volkoff limit (maximum mass of neutron star) M ^ 1.5 or 2 M 

ov e 

To verify for a given star that giant-type evolution (unchanging 
envelope abundances) really is reasonable, and to learn the details of 
evolution of the star's deep interior, one must construct a sequence of 
evolutionary models for the entire stellar envelope and halo. One could 
do so using a Henyey-type code. However, we think this is not necessary. 
Assuming that our models are stable, Henyey evolution would have to re- 
produce the unique sequence of models which we obtain by our methods holding 
the star's total rest mass M fixed and increasing its core mass from 
model to model. Such a sequence would be nearly the same as the one shown 

in Table 4 with fixed total mass M and increasing M. . 

t tc 

The evolutionary sequence in Table 4 is for a star with total mass 

M = 5 M , with envelope abundances X, =0.70, Y , =0.27, C. =0, and 
to’ r ph ph ph 

with core mass increasing from 0.40 M initially to a final, Oppenheimer- 

0 

Volkoff limit of 1.625 M . The core radius-mass relation R (M ) is that 

© c tc 

of Malone, Johnson, and Bethe (1975, their model V-H) . In our models we 
were satisfied with reproducing the desired to within about one 

percent. All models in our evolutionary sequence (Table 4) have giant-type 
structures. As one might expect, as the core mass grows and the accelera- 
tion of gravity at the core edge increases, the thickness of the halo de- 
creases (cf. §IV.b). The total time required for evolution from M = 

0.4 M to the point of core collapse, M = 1.625 M , is 
o ‘ tc 0 
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(7.1) 


At - M ^dM = -7.4 XIO^ years . 

J c J c tc J 

The evolution of a supergiant is more complex than that of a giant; 
it is driven not only by core growth, but also by chemical evolution of 
the envelope. The rate at which envelope hydrogen is burned by the shell 
sources of our models to form envelope helium is typically ^500 times 
greater than the rate at which envelope matter flows into the core; see 
Table 3. To burn all of its envelope hydrogen a supergiant of 12 re- 
quires ^1.1x10^ years, and a supergiant of 25 M requires VL.4x 10^ years. 

0 

For comparison, the time required for the core rest mass to increase 

8 8 

1 M # is ^ 6 x 10 years in the first case and ^7x10 years in 

the second. These estimates may be in serious error because they are based 

on our inadequate treatment of the hydrogen burning (§V) and on models 

(Table 3) of one chemical composition only. We have not yet attempted to 

construct supergiant models with hydrogen-deficient envelopes. 

VIII. CONCLUSION 

We regard this paper as merely a first rough overview of stellar models 
with neutron cores. This ove ~w has uncovered a large number of problems 
which must be resolved before the theory will be in satisfactory shape. At 
present we are pursuing only three of these problems vigorously : the details 
of nuclear burning and nucleosynthesis in supergiant models (§V) , the re- 
sulting chemical evolution (§VII), and the possibility of discovering such 
stars by observation of peculiar photospheric abundances (§I.c), 

Other problems that require study are these: (i) The stability of our 
models, with emphasis on the five possible instabilities described in §VI. 
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(ii) A search for models with very different structures from those ex- 
hibited in this paper — e.g., models with large neutrino losses supplied 
by large mass inflow rates (cf.§IV.a and §VI) and models in which pro- 
fuse electron-positron pairs replace large as the source of reason- 
able scale heights above the halo (§IV.a). (iii) The effect of convec- 
tive overshoot on the structures of our giant models (§IV.b). (iv) The 

effect of magnetic fields in reducing the opacity at densities 
6 11 

10 < p < 3 x 10 and thereby permitting significant heat transfer between 

core and halo (§11. b). 

All of the above problems seem somewhat tractable. Less tractable, 
but obviously very important, is the issue of how such stars might form 
in Nature (Ostriker and Paczynski 1975). 
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APPENDIX 


ELECTRON-POSITRON PAIRS IN THE NONRELATIVISTIC, NONDEGENERATE APPROXIMATION 

In thermodynamic equilibrium at temperature T , the number density of 
electrons and positrons "+" in phase space is 


n- = 


dNq. 


+ d 3 xd 3 p 


h 3 1 + exp[(p Q + y)/kT] 


(A. 1) 


2 2 1/2 

where y is the chemical potential, p^ = (m + p ) is the total mass- 
energy of a particle, p is the magnitude of the spatial part of its 
4-momentum, m is the electron rest mass, k is Boltzmann’s constant, 
and we set the speed of light equal to unity. Here and below equations con- 
taining double signs (+ or +) are valid for electrons (including ionization 
electrons and pair electrons) with the upper sign and for positrons with the 

lower sign. The number densities in physical space n , the pressures P , and 

* + 

the energy densities including rest mass are 


n f .]V f * 3 P , P— i 


(p 2 /p 0 )^ _ d 3 p , e_ = 
u + + 


P 0 ^ d P > (A- 2) 


3 2 

where d p = 4TTp dp . 


We assume that the number of electrons exceeds the number of positrons, 
so y > 0 ; and we specialize to the nonrelativistic, nondegenerate regime: 


y > 0 


m/kT » 1 


(m-y ) /kT » 1 


(A. 3) 


In this regime relativistic particles make negligible contributions to 

n_, P_, and j[ ; consequently, we can set p /2m << m , and use 

+ + + 

2 

P 0 = m + p /2m . (A. 4) 

Assumptions (A. 3) then allow us to write 
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+ h" 


71. - \ exp(- e X p( i,J -- ra - 


kT 


) 


(A.5) 


By inserting expressions (A. 4) and (A.5) into (A. 2) and integrating 
we obtain the following results: (i) The electron and positron pressures 
and energy densities are given by the usual nonrelativis tic expressions 


P_ = n_ kT , - n_(m + — kT) 

+ + + + 2 


(A. 6a) 


(ii) The number densities of electrons and positrons are 


n_ = (2/h 3 ) (2T[mkT) 372 exp[ (±y - m) /kT] 


(A. 6b) 


mean 


The number density of ionization electrons n g is the difference 
between n and n + > — and is also equal to p/m^p^ where p g is the 
molecular weight per electron: 


p/m H P e = n g = n_ -n + = (2/h 3 ) (2TrmkT) 3/2 e m/kT [ e y/lcT - e p/kT ].(A.7) 


Let us introduce the parameter 

m/kT / P 6 /U e\/5.93\ 3 / 2 5.93/T, 

" \ 7. 37/ \ T 


i 3 
h n 


y = 


4(27imkT) 


3/2 




9 I 


(A. 8) 


6 3 

where Pg is density in units of 10 g/cm and T^ is temperature in units 

9 li/kT 

of 10 K. Then by solving equation (A. 7) for e M/ we obtain 


P/kT , , 2, .,1/2 
e = y + (y + 1) 


(A. 9) 


and by combining with (A. 6b) we obtain for the ratio of number of electron- 
positron pairs to number of ionization electrons 


n . 


n 


n - n 


n 


2y[y + (1+y 2 ) 172 ] 


(A. 10) 
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Equation (A. 6a) then shows that the ionization electrons make the usual 
contribution to the pressure and energy density, while the pairs make the 
contribution (in cgs units) 


P . 
pair 


pn 


pair 


2n + kT 


(P/y^kT 

y[y +(i+y 2 ) 1/2 l 


2n + (mc^ + —■ kT) 


(p/y e m H )(mc2+ J 
y [y + (l+y 2 ) 1 ^ 2 ] 


(A.lla) 


(A. lib) 


The temperature-density regime in which the above expressions are 
valid can be deduced by combining equations (A, 8) , (A. 9), and (A. 3): 


T 9 « 5.93 


3.69 y 


« 


T 9 \ 3/2 
5 .93 / 


(A. 12) 


The T-p curve along which the number of pairs equals the number of ioni- 
zation electrons is given by y = 0.354 and is shown graphically in Figure 
2. Note that when our stellar interiors cross over this curve so pairs 
become important, they remain well within the realm of validity of our 
approximations, equations (A. 12)* 
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TABLE 1 


M tc = 1 -° V R c = 10 - 00 X ph = °- 70 > Y ph = °' 27 ' C ph = °*° 

. r . = 9.95 km. For details of notation see §I.e.) 
r m-i 


T 

L 

r 

(V 

l ov -S 2 l ov 

r 

(V 

L cr it 

r 

(V 

V ad 

V , 
rad 

V 

log v t 
(cra/sec) 

iog v s 
(cm/ sec) 

l0g V in 
(cm/sec) 

0.667 




0.0834 

0.106 

0.106 

2.34 

5.63 


17.9 




.386 

1.63 

1.63 

4.43 

5.87 


8970. 




.0869 

168. 

0.519 

5.42 

5.92 






.210 

298. 

.534 

5.54 

6.25 






.283 

6.57 

.496 

5.89 

6.60 






.2502 

0.340 

.264 

6.60 

7.31 






.2500 

.2775 

.2516 

7.18 

7.92 


l-fi^/L 










0 

41745. 

« E-20 

38160. 

.2500 

.275 

.2504 

7.64 

8.40 

5.00 

-0.010 

42306. 

« E-20 

38640. 

.2500 

.274 

.2500 

8.12 

9.02 

6.46 

- .034 

44127. 

« E-20 

40008. 

.2500 

.275 

.2500 

8.27 

9.35 

6.85 

- .058 

47348. 

2. IE-20 . 

43061. 

.2500 

.274 

.2500 

8.25 

9.56 

6.81 

- .073 

52695. 

3.5E-18 

49306. 

.2500 

.266 

.2500 

8.12 

9.71 

6.60 

- .078 

61858. 

4.0E-16 

61246. 

.2500 

.2522 

.2500 

7.75 

9.84 

6.32 

- .078 

63501. 

7.4E-16 

63449. 

.2500 

.2503 

.2500 

7.44 

9.85 

6.28 

- .078 

63662. 

7.8E-16 

63666. 

.2500 

.2501 

0.2500 

7.28 

9.85 

6.28 







X 

Y 

C 


- .078 

63741. 

8.0E-16 

63772. 

.2500 

.2500 

0.700 

0.270 

0 

6.27 

- .039 

61440. 

8.2E-16 

63829. 

.2500 

.2423 

.700 

.270 

0 

6.26 

+ .456 

32201. 

8.2E-16 

63845. 

.250 

.138 

.700 

.270 

0 

5.96 

.709 

17201. 

8.2E-16 

63850. 

.250 

.0764 

.700 

.270 

0 

5.67 

.913 

5126. 

8.4E-16 

63860. 

.250 

.0227 

.700 

.270 

0 

4.99 

.966 

2025. 

1.2E-15 

63897. 

.251 

.0095 

.700 

.270 

0 

3.72 

.969 

1836. 

1.1E-14 

69345. 

.284 

.0296 

.700 

.270 

0 

2.45 

.970 

1748. 

2.9E-13 

69090 

.348 

.129 

.669 

.301 

0 

1.43 

.982 

1049. 

1.5E-12 

90198. 

.358 

.094 

.328 

.642 

0 

0.94 

.9927 

434. 

3.0E-12 

1.18E5 

.362 

.041 

0.026 

.944 

0 

0.59 

.9936 

377. 

1.2E-11 

1.23E5 

.371 

.075 

0 

.970 

0 

0.19 

.9938 

365. 

9. IE- 11 

1.33E5 

.372 

.141 

0 

.930 

0.040 

-0.29 

.9955 

265. 

4.9E-10 

1.48E5 

.363 

.135 

0 

.441 

.529 

-0.62 

.9969 

181. 

8.3E-9 

2.01E5 

.343 

.139 

0 

0.030 

.940 

-1.11 

.9971 

172. 

8.6E-7 

4.07E5 

.322 

. 158 

0 

0 

.970 

-1.79 

.9974 

156. 

8.5E-5 

9.17E5 

.312 

.141 

0 

0 

.894 

-2.44 

.9988 

74. 

2.5E-4 

1.19E6 

.313 

.069 

0 

0 

.407 

-2.60 

0.9998 

9. 

1.2E-3 

2.07E6 

.318 

.011 

0 

0 

.035 

-2.87 

1.0000 

3. 

4.6E-3 

3.73E6 

.325 

.004 

0 

0 

0.005 

-3.12 

1.0000 

0. 

1.9E-2 

7.49E6 

0.336 

0.000 

0 

0 

0 

-3.41 



TABLE 2 





INTERNAL STRUCTURE OF 

A STAR 

WITH M = 

12-0 V 

M. a 1.0 
tc 

II 

sP 

9.98 km, : 

X , = 0.7C 
ph 

Y P h = 

: 0.27, C 

, = 0.0 
ph 







(Note: r^, = 

= 10.00 

tab M tR = 

1.000 

r . = 9. 

m-i 

.92 kmj for details 

of notation see 

§I.e.) 






Region 

r ' r K 

M tr"V 

2GM tr 

2 

log p 

log T 

p g 

T 

L 

r 

l ov -R 2 l° v 

L crit 

r 

7 ad 

^rad 

V 

log v t 

iog v s 

106 v in 



(m) 

(H,) 

C T 

(g/cra 3 ) 

(°K) 



(V 

<v 




(cm/ sec) 

(cm/sec) 

(cm/sec) 

PHOTOSPHERE 

7.50E11 

10.90 

4.69E-8 

-8.397 

3.462 


0.667 




0.118 

0.148 

0.148 

2.69 

5.66 


i 


7.43E11 

10.80 

4.69E-B 

-8.191 

3.598 


3.92 




.387 

1.01 

1.01 

3.66 

5.81 




7.3BE11 

10.68 

4.68E-8 

-8.330 

3.916 


1590. 




.104 

76.5 

0.836 

5.45 

5.90 




6.71E11 

9.44 

4.60E-8 

-8.127 

4.200 






.232 

470. 

.629 

5.53 

6.21 




4 -03E11 

5.58 

3.36E-8 

-7.709 

4.671 






.254 

23.9 

.352 

5.55 

6.49 


' STATIC 

3.28E10 

S.45E-3 

9.06E-8 

-7.194 

5.148 






.263 

2.55 

.352 

6.09 

6.05 


i ENVELOPE , 

2.15E9 

5.06E-6 

1.37E-6 

-6.699 

5.617 






.250 

0.594 

.265 

6.64 

7.38 




1.35E8 

6.84E-9 

2.19E-5 

-5.999 

6.092 






.250 

• 501 

.252 

7.19 

7.95 










l-R 2 L r /L 










X 

r 


1.53E7 

6.90E-11 

1.93E-4 

-5.243 

6.504 

9.21E-3 

0 

73547. 

« E-20 

38301. 

.250 

.484 

.251 

7.57 

8.40 

2.79 

c 


1.48E6 

2.33E-12 

1.98E-3 

-3.996 

7.040 

4.04E-3 

0 

73691. 

6.1E-21 

39000. 

.250 

.474 

.250 

7.82 

8.84 

3.57 


i . 

4.04E5 

9.70E-13 

7.12E-3 

-2.963 

7.417 

3.23E-3 

0 

74086. 

9.1E-19 

40411. 

.250 

.459 

.250 

7.85 

9.08 

3.65 



1.45E5 

6.22E-13 

0.0191 

-1.958 

7.762 

3.01E-3 

-0.001 

75003. 

2.0E-16 

43448. 

.250 

.431 

.250 

7.78 

9.26 

3.50 

1 INFLOWING 1 

5.64e4 

4.32E-13 

.0445 

-0.969 

8.095 

2.95E-3 

-0.001 

77003. 

4.8E-14 

50043. 

.250 

.382 

.250 

7.67 

9.43 

3.25 

, ENVELOPE , 

2.42E4 

3.01E-13 

.0864 

-0.119 

8.379 

2.93E-3 

-0.001 

80545. 

5.8E-12 

61612. 

.250 

.324 

.250 

7.53 

9.57 

2.99 



1.03E4 

1. 93 E- 13 

.146 

+0.596 

8.617 

2.93E-3 

0 

86042. 

4.0E-10 

79400. 

.250 

.270 

.250 

7.28 

9.68 

2.74 

1 

H 

Shell begins 

4 .4 5E3 

1.12E-13 

.204 

1.097 

8.784 

3.96E-3 

+0.024 

90228. 

6.4E-7 

59575. 

.249 

.366 

.249 

7.42 

9.76 

2.55 



2.51E3 

7.29E-14 

.236 

1.324 

8.858 

8.84E-3 

.100 

86617. 

1.5E-5 

23150. 

.247 

0.857 

.247 

7.48 

9.78 

2.46 



883. 

2.97E-14 

.271 

1.559 

8.932 

0.0217 

.401 

60407. 

2.0E-4 

9740. 

.243 

1.396 

.243 

7.40 

9.81 

2.36 



126. 

4.61E-15 

.292 

1.688 

8.972 

.0334 

.831 

17518. 

6.7E-4 

6719. 

.241 

0.612 

.241 

7.16 

9.82 

2.30 



14.4 

5.34E-16 

.295 

1.708 

8.9782 

.0356 

.931 

7222. 

8.0E-4 

6377. 

.240 

.288 

.240 

6.85 

9.82 

2.29 



5.56 

2.07E-16 

.295 

1.710 

8.9787 

.0357 

.939 

6347. 

8.1E-4 

6351. 

.240 

.259 

0.240 

X 

6.71 

Y 

9.82 

c 

2.29 

V 

K1 

i 

JEE 

0 

0 

.295 

1.711 

8.9790 

.0359 

.945 

5788. 

8.2E-4 

6336. 

.240 

.240 

0.700 

0.270 

0 

2.28 

/ 

f 

-0.57 

-2.54E-17 

.295 

1.905 

8.9791 

.0360 

.964 

3705. 

8.2E-4 

9817. 

.240 

.114 

0.367 

0.603 

0 

2.09 

H 

Shell ends 

-1.13 

-7.36E-17 

.295 

2.302 

8.9791 

.0365 

.986 

1485. 

8.2E-4 

24194. 

.240 

.0267 

0.021 

0.949 

0 

1.69 



-3.82 

-1.06E-15 

.296 

3.045 

8.9791 

.0491 

.9971 

301. 

8.3E-4 

1.16E5 

.241 

.0029 

0 

0.970 

0 

0.95 

lie Shell i 

( -8.59 

-1.07E-14 

.296 

3.819 

8.9792 

.159 

.9992 

85. S 

8.3E-4 

2.19E5 

.247 

.0005 

0 

0.949 

0.021 

0.18 


1 ■ ■ . 

-12.5 

-4.62E-14 

.296 

4.453 

8.9792 

.411 

.9996 

42.5 

8.4E-4 

2.28E5 

.253 

.0002 

0 

0.480 

0.490 

-0.46 

' HALO 1 1 

[-17.6 

-3.70E-13 

.296 

5.385 

8.9793 

.843 

.9998 

22.1 

8.5E-4 

2.29E5 

.283 

.0002 

0 

0.022 

0.94a 

-1.39 


f 

-19.7 

-8.66E-15 

.296 

5.732 

8.9794 

.926 

.9998 

20.1 

8.5E-4 

2.30E5 

.299 

.ooo4 

0 

0.006 

0.959 

-1.73 


1 

[ -22.1 

-2.08E 12 

.296 

6.068 

8.9796 

.968 

.9998 

18.4 

8.6E-4 

2.36E5 

.309 

.0007 

0 

0.002 

0.940 

-2.07 

c 

Shell j 

-27. B 

-1. 15E-11 

.296 

6.667 

B.9804 

.994 

.9999 

9.5 

8.8E-4 

2.83E5 

.310 

.0015 

0 

0 

0.491 

-2.67 


! 

‘ -33.2 

-6.33E-11 

.296 

7.220 

8.9812 

0.999 

1.0000 

1.2 

9.7E-4 

5.01E5 

.306 

.0007 

0 

0 

0.032 

-3.22 

‘‘ core ‘ 

' -49.7 

-3.S0E-10 

0.297 

7.776 

8.9816 

1.000 

1.0000 

0.1 

1.9E-3 

1.44E6 

0.312 

0.0003 

0 

0 

0 

-3.78 


TABLE 5 


MODELS WITH M t( . = 11^, R c = 10.00 km, X ph =0.70, Y ph =0.27, C ph = 0.0 


Type M 

Model* L 

(Mj 

L 

(V 

1 

ph 

(°K) 

R ph 

108 Vm 

(°k) 

M 

(10' 8 1^/yr) 

L /L 
nuc' 

r K 

(km) 

108 P K 
(g/ cm 3 ) 

log T k 
(°K) 

r H -I k 

(m) 

log P H 
(g/cm 3 ) 

iog T„ 
(°K) 

r He- r K 

(m) 

108 p He 
(g/cm 3 ) 

108 T He 
(°K) 

r C~ r K 
( m ) 

log P c 

(g/cm -5 ) 

Envelope 

Stable?* 

G,R,0: S3 1 2?.»Q 

38220 

2621 

917 

5.923 

1.575 

0.030 

10.039 

-1.871 

8.098 

-37 

4.12 

8.369 

-4l 

5.61 

8.529 

-66 

7.62 

NO 

G,R,Q=1 3.0 

38430 

2549 

984 

5.977 

1.584 

0.030 

10.045 

-1.817 

8.111 

-36 

4.12 

8.369 

-40 

5.61 

8.530 

-65 

7.62 

NO 

G,R,Or=l 5.0 

41740 

2649 

961 

6.190 

1.720 

0.030 

10.047 

-1.267 

8.249 

-37 

4.05 

8.386 

-4l 

5.60 

8.537 

-65 

7.59 

YES 

G,R.,a=l 8.0 

55680 

2776 

1016 

6.374 

2.286 

0.030 

10.024 

-0.004 

8.561 

-47 

3.35 

8.569 

-53 

5.41 

8.610 

-75 

7.46 

YES 

G,R,Ct = 1 9.0 

MASS 
GAP 

60380 

2809 

1034 

6,406 

2.477 

0.030 

10.028 

0.309 

8.639 

-58 

3.10 

8.643 

-65 

5.29 

8.663 

-87 

7.41 

YES 

S,R,a«l 11.5 

71360 

2882 

1070 

6.498 

0.167 

0.945 

10.000 

1.706 

8.978 

711 

1.58 

8.940 

-12 

4.44 

8.978 

-28 

6.67 

YES 

S,R,a = l 12.0 

73530 

2894 

1077 

6.504 

0.165 

0.947 

10.000 

1.711 

8.979 

725 

1.59 

8.940 

-12 

4.44 

8.979 

-28 

6.66 

YES 

S,R,a=l 16.0 

90490 

2978 

1130 

6.539 

0.156 

0.959 

10.000 

1.743 

8.987 

712 

1.62 

8.949 

-12 

4.42 

8.987 

-27 

6.57 

YES 

S,R,a=l 20.0 

105600 

3047 

1168 

6.563 

0.150 

0.966 

10.000 

1.766 

8.992 

731 

1.64 

8.953 

-12 

4.4l 

8.992 

-26 

6.51 

YES 

S,R,Q = 1 25.0 

123260 

3114 

1209 

6.585 

0.144 

0.972 

10.000 

1.769 

8.998 

718 

1.66 

3.960 

-11 

4.4o 

8.998 

-25 

6.45 

YES 

G,N,Q - 1 5.0 

41770 

2649 

961 

6.179 

1.852 

0.038 

10.040 

-2.071 

8.026 

-32 

4.07 

8.369 

-38. 

5.60 

8.532 

-72 

7.59 

YES 

G,R,a = 0.5 5.0 

R,Q=1.5 5.0 

38180 
IN MASS 

2209 

GAP; 

1324 6.029 

NO MODEL EXISTS 

1.573 

0.030 

10.068 

-1.693 

8.145 

-62 

4.11 

8.370 

-66. 

5.63 

8.526 

-93 

7.69 

YES 

* 

In the first column, G means giant. 

S means supergiant, R means relativistic model, N 

means Newtonian mode, and CC is 1 

Che ratio 

of mixing length 

to pressure scale height. 

In the last column "ENVELOPE STABLE? 

" means "Is the static envelope stable against small, adiabatic, 

TABLE 4* 

EVOLUTIONARY SEQUENCE FOR A STAR OF M fc = 5}^, X ph = 0. 

radial perturbations?" For other 
70, Y ph = 0.27, C ph = 0.0 

details of notation see 

§I.e. 

M R 

tc c 

L 

T t, 
ph 

R u 

ph 

log T „ m 

o-m 

M 

L /L 
nuc' 

r K 

iog p K 

!ogT K 

V r K 

!og P H 

108 t h 

r He" r K 

iog P Hg 

log T He 

r C -r K 

log P c 

Envelope 

(H.) (km) 

(V 

( q k) 

(id) 

(°K) 

do' 8 n,/y r ) 


(km) 

(g/ cm 3 ) 

(°K) 

(m) 

(g/cm 3 ) 

(°K) 

(m) 

(g/cm 3 ) 

(°K) 

(m) 

(g / cm 3 ) 

Stable?* 

0.1 12.172 

24590 

2734 

696 

6.569 

3.028 

0.100 

12.398 

0.421 

8.524 

-181 

3.07 

8.54 

-209 

5.16 

8.62 

-273 

6.94 

YES 

0.8 11.255 

36210 

2659 

889 

6.298 

2.121 

0.045 

11.323 

-0.975 

8.275 

-57 

3.87 

8.40 

-64 

5.48 

8.55 

-98 

7.36 

YES 

1.0 11.075 

41730 

2649 

961 

6.193 

1.913 

0.034 

11. 124 

-1.529 

8.219 

-41 

4.01 

8.38 

-46 

5.56 

8.54 

-75 

7.52 

YES 

1.1 10.680 

54400 

2654 

1091 

5.999 

1.673 

0.021 

10.710 

-1-505 

8.229 

-26 

4.16 

8.37 

-29 

5.68 

8.52 

-50 

7.82 

YES 

1.6 10.250 

61860 

2670 

1140 

5.956 

1.561 

0.016 

10.277 

-1.329 

8.300 

-24 

4.16 

8.38 

-26 

5.73 

8.52 

-43 

7.92 

YES 

1.625 9.299 

62810 

2672 

1156 

5.950 

1.388 

0.013 

9.336 

-1.153 

8.364 

-35 

4.10 

8.41 

-37 

5.75 

8.53 

-49 

7.93 

YES 


See footnote Co Table 5. 
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FIGURE CAPTIONS 


Figure 1. The structure of stars with degenerate neutron cores. The 

interior of the box lists a number of features of the stellar interior. 
The locations of those features are indicated on the left of the box 
in terms of density p , and on the right of the box in terms of radius 
minus the radius of the star’s "knee", r-r^ ♦ Numbers not in parenthe- 
sis are exact and apply to all of our models. Numbers in parenthesis are 
taken froisi a general relativistic "giant" model (Table 1) with total mass 

M = 5 M , total core mass M = 1 M , and core radius R = 10 km — but 
to tc ® u 

these parenthesized numbers are qualitatively correct for all models. 

In the left column of the box are listed the major regions of 
the stellar interior: the photosphere , which is the point with optical 
depth T = 2/3; the envelope , which extends downward from the photosphere 
to the knee; the knee , which is the point where envelope convection 
stops; the halo , which extends inward from the knee to the point of 
onset of electron degeneracy; and the core , which extends from the onset 
of electron degeneracy in to the center of the star. 

In the middle column of the box are listed a number of subregions 
of the stellar interior including: the atmosphere , which lies above 
the photosphere; the static envelope , which is a convective region ex- 
tending from the photosphere down to a (arbitrarily chosen) radius 

r where inflow of matter from envelope to core becomes important; 
o-m 

the inflowing envelope , which is also convective and extends downward 

from r to the star's knee where convection ceases; the region of 
o— m 

gravitational-energy release , which extends inward from the knee to 
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a density p ^ 10p^ ; ^Y^ro gen-burning shell , helium-burning 

shell , and carbon-burning shell ; an insulating layer which extends 

from the onset of electron degeneracy down to the point p = 

11 3 

3x10 g/cm where neutrons drip off the atomic nuclei to form a 

superconducting, superfluid medium; and the isothermal core in which 
X / 2 

T I ^ 00 1 ~ ^ i- s nearl y constant, and which extends from 

11 

p - 3 x 10 in to the center of the star. 

In the right column are listed three regions into which we sub- 
divide the model for computational purposes: The outer region , which 
includes atmosphere, photosphere, and static envelope; the middle region , 

which includes contracting envelope, halo, and the outermost part of 

8 3 

the core ("down to radius r . where D - 3x10 g/cm ); and the inner 

m-i 

region , which includes the remainder of the core. 

Supergiant models (M > 10 M ) differ from giant models (M < lOM^; de- 
picted in this figure) in only one qualitative way; the hydrogen-burning 
shell overlaps the knee and envelope instead of being confined to the 
halo. 

Figure 2. The internal distributions of density and temperature for a 

giant model with total mass M = 5 M , and a supergiant with M = 
a to t 

12 M . Both models have envelope abundances (= photospheric abund- 
0 

ances) X , = 0.70, Y , = 0.27, C , = 0, and core mass and radius 
ph ph ph 

M = 1 M , R = 10 km. Further details of the internal structures of 
tc 0 c 

these stars are given in Tables 1 and 2. The solid curves are the runs 
of density and temperature in the two models . Along these curves are 
marked several regions of the model which were described qualitatively 
in Figure 1 (photosphere, static envelope, junction point r m between 
outer and middle regions, inflowing envelope, knee, halo, junction 
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point R^ between halo and core, core, and junction point r^ ^ 

between middle and inner regions. Also shown along each curve are 

the locations of the hydrogen-, helium-, and carbon-burning shells. 

The dotted lines are several regions of interest in the density- 

temperature plane: Above the "PAIRS — NO PAIRS" line the density of 

electron-positron pairs exceeds that of ionization electrons; below, 

ionization electrons dominate. (This curve is given analytically by 

y = 0.354 where y is defined by eq.(A.8) of the Appendix). The 

"RADIATION DOMINANCE-- GAS DOMINANCE" line is the line where g = 

g 

(gas pressure) /(radiation pressure plus gas pressure) is 0.5. For 
further details on notation see §I.e. 

Figure 3. H-R diagram (luminosity versus photospheric temperature) for 

stars of total mass M = 5 M , envelope abundances X , = 0.70, Y = 0.27, and 

t 0 ph ph 

ratio of mixing length to pressure scale height ^ t / H p = 1* Each point in the 
L-T ^ diagram corresponds to a unique static envelope constructed by the pre- 
scription of §II.c. Each dotted curve is a region of constant core mass, 

M - M . The curve M = 0 was calculated by extrapolation from M^ > 0. Each 
t env tc tc 

thick solid curve is a region of constant temperature T Q _ m at the 
base of the envelope (radius r 0 _ m defined by eq. [2.1]), To the 
left of the thin solid curve the turbulent velocity of convection is 
less than half the adiabatic sound velocity throughout the static 
envelope. To the right, v £ > v g /2 near the base of the envelope — 
and therefore we are not justified in our use of subsonic mixing-length 
theory. The large dots are the surface properties of the six stellar 
models of Table 4. 

The absolute temperatures and luminosities along the various 
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curves are unreliable because of uncertainties in opacities and 
mixing length (which is here assumed equal to pressure scale 
height) . However, the temperature and luminosity differences be- 
tween various curves should be somewhat reliable. 

Figure 4. Local luminosity and critical luminosity L^ rit as func- 

tions of temperature T in the inflowing envelopes of the 5 

giant model of Table 1, and the 12 M q supergiant model of Table 2. 

ciri t 

We actually plot vertically and multiplied by the red- 

G rit 

shift factor R because the product RL^ is very nearly a 
function of temperature only and is therefore the same for all 
models with the same core masses (cf. eqs . [5 . 1] , [5. 2] ) . The L^ 
curves are parametrized by radius r in kilometers. The knee 
(r = 10.0 km) occurs where goes subcritical. 
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